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SECTION  1 
INTRODUCTION 

It  is  well  known  that  the  mechanical  properties  of 
materials  may  be  influenced  by  strain  rate.  Many  materials  do 

•  not  exhibit  strain-rate  dependency  at  quasi-static  strain-rate 
levels.  However,  most  structural  materials  show  rate  sensitivity 

*  above  a  particular  strain-rate  level,  typically  100/s.  The 
magnitude  of  strain-rate  effects,  especially  in  metals,  has  been 
a  topic  of  numerous  investigations. 

Several  applications  of  structural  materials  involve 
impulsive  loading.  Resulting  deformations  are  often  very 
complicated.  The  material  behavior  may  involve  a  high  rate  of 
strain,  large  deformation,  high  pressure  and  high  temperature  due 
to  adiabatic  heating.  Dynamic  materials  characterization  is 
especially  important  in  the  analysis  of  weapons  effects;  problems 
include  structural  response  to  explosive  loading  as  well  as  local 
failure  from  impact  of  gun-launched  and  explosively-launched 
projectiles . 

Within  the  last  decade,  the  development  of  finite 
element/difference  codes  has  provided  additional  analytical 
capability  in  understanding  these  problems.  It  is  now  well 
established  that  the  material  descriptions  can  affect  the 
computational  results  very  significantly  (see,  for  example 
Reference  1).  In  most  computations,  the  material  behavior 
parameters  have  been  indirectly  determined  by  adjusting  input 
parameters  in  order  to  obtain  agreement  with  experimental 
‘  observations  (a  process  known  as  "post  shot  prediction").  This 

indirect  method  of  obtaining  material  properties  such  as  dynamic 
'  yield  strength  can  often  be  misleading  since  error  in  the 

physical  model  can  be  masked  by  unrealistic  material 
descriptions.  Since  dynamic  material  property  data  are  not 
readily  available,  and  there  are  no  simple  tests  from  which  to 
obtain  these  data,  code  users  have  often  been  unable  to  apply 
sophisticated  constitutive  relations  for  describing  materials. 
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For  example,  most  computer  codes  still  do  not  account  for  the 
strain-rate  and  pressure  dependence  of  yield  and  flow  stresses. 

Nicholas  (Reference  2)  describes  in  detail  the  various 
experimental  techniques  that  are  being  currently  used  by  several 
investigators  to  characterize  materials  under  dynamic  loading 
conditions.  The  present  report  describes  the  combined 
experimental  and  theoretical  efforts  undertaken  to  model  the  high 
strain-rate  material  response  of  1020,  HYlOO,  C1008  steels,  OFHC 
copper,  7039-T64  aluminum,  and  BeO  ceramic.  This  work  is  an 
extension  of  the  results  presented  by  Rajendran  et  al., 

(Reference  3).  The  report  also  describes  the  development  of  new 
techniques  for  high  strain  characterization:  a  double  flyer 
plate  experiment  to  study  spall  and  recompaction,  and  a  Cranz- 
Schardin  camera  to  observe  transient  profiles  in  Taylor  impact 
specimens . 

Dynamic  tensile  and  compressive  loading  under  a  one¬ 
dimensional  stress  state  was  achieved  with  a  split  Hopkinson  bar 
(SHB).  SHB  data  were  extended  to  higher  strains,  strain-rates, 
and  mean  stresses  by  high  speed  photography  of  necking  specimens. 
The  plate-impact  test  provided  spall  threshold  data  and  yield 
strength  at  very  high  strain-rate  levels. 

The  state  variable  based  visco-plastic  constitutive  theory 
of  Bodner  and  Partom  (Reference  4)  was  used  to  characterize  the 
metals  investigated  in  this  program.  A  series  of  automated 
computer  programs  were  developed  to  evaluate  the  model  parameters 
from  SHB  (Tension)  and  plate  impact  test  results.  The  SHB 
tension  and  compression  tests  were  simulated  with  the  developed 
model  parameters  and  compared  with  the  actual  test  results.  The 
plate  impact  tests  were  numerically  simulated  through  a  state-of- 
the-art  general  purpose  finite  difference  code,  called  'STEALTH', 
using  the  Bodner-Partom  (BP)  model  parameters  for  each  material. 
For  this  purpose,  special  purpose  subroutines  describing  the  BP- 
model  were  developed  for  STEALTH. 

Two  failure  criteria  were  considered.  The  first  criterion 
was  time  independent  and  based  on  a  critical  spall  stress.  The 
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spconc^  was  based  on  a  critical  value  for  a  time  dependent 
integral  as  proposed  by  Tuler  and  Butcher.  A  zone  model  for 
spall  in  ceramics  was  also  developed  and  verified.  Failure  model 
parameters  were  obtained  from  the  numerical  simulations  of  the 
plate  impact  tests. 
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SECTION  2 

EXPERIMENTAL  PROGRAM 

The  experimental  techniques  considered  for  investigating  high 
strain  rate  material  response  were  the  split  Hopkinson  bar  and  the 
one-dimensional  plate  impact  test  configurations.  These  two  tests  » 

encompass  extremes  of  stress  and  strain  states.  The  conventional 
SHB  test  exerts  a  uniaxial  stress  state  in  which  all  three  strain  * 

components  are  non-zero.  On  the  other  hand,  plate  impact  leads  to 
a  uniaxial  strain  state  in  which  all  three  stress  components  are 
non-zero.  The  experimental  program  also  included  an  unconventional 
SHB  test  which  employed  high  speed  photography  of  necking 
instabilities.  This  'necking  SHB  test'  provides  an  effective 
method  to  extend  the  SHB  data  to  higher  strains,  strain-rates,  and 
mean  stress.  In  this  section,  the  experimental  techniques  are 
described  and  the  principal  results  are  summarized. 

A  parallel  activity  on  this  contract  has  been  development  of 
new  experimental  techniques  for  high  strain  rate  characterization. 

Two  approaches  have  been  pursued.  The  first  is  a  double  flyer 
plate  technique  to  study  recompaction  of  voids.  The  second  is  use 
of  a  Cr anz-Scha rd in  camera  to  resolve  intermediate  contours  in 
Taylor  impact  tests.  Results  of  these  development  efforts  are 
reported  in  Section  2.4. 

2.1  SPLIT  HOPKINSON  BAR 

The  split  Hopkinson  bar  provides  one  of  the  few  research 
tools  for  investigating  the  behavior  of  materials  under  uniaxial 
stress  loading  at  strain  rates  above  300  s  The  University's 

Hopkinson  bar  has  been  designed  to  measure  tensile,  as  well  as 
compressive,  stress-strain  relationships. 

2.1.1  Technique 

A  split  Hopkinson  bar  consists  of  three  in-line  bars, 
a  striker  bar,  a  pressure  bar,  and  a  transmission  bar,  as  shown  in 
Figure  1.  The  University's  bars  are  12.7mm  in  diameter.  The 
striker  bar  is  launched  by  a  torsional  spring.  Its  speed  is 


4 


Figure  1. 


Hopkinson  Bar  Apparatus. 


measured  as  it  crosses  two  lamp/photodetector  stations.  The 
striker  bar  strikes  the  pressure  bar,  producing  an  elastic  com¬ 
pressive  stress  wave  traveling  at  Cg=  / E/p ,  where  E  is  Young's 
modulus  and  p  the  mass  density  of  the  bar.  The  duration  of  the 
stress  wave  is  twice  the  transit  time  through  the  striker  bar, 

300  ps.  For  compressive  tests,  a  button  sample  is  placed  between 
the  pressure  and  transmission  bars.  For  tensile  tests,  a  collar  is 
placed  around  the  specimen  to  transmit  the  compressive  pulse  from 
the  pressure  bar  to  the  transmission  bar.  The  pressure  wave 
reflects  at  the  free  end  of  the  transmission  bar  and  returns  as  a 
tensile  wave.  The  specimen,  shown  in  Figure  2,  is  screwed  into 
both  the  transmission  bar  and  the  pressure  bar.  The  specimen  is 
loaded  in  tension  by  the  tensile  pulse  arriving  in  the  transmission 
bar.  Analyses  of  the  stress  waves  reflected  and  transmitted  by  the 
sample  can  be  used  to  deduce  the  stress-strain  history  of  the 
specimen.  The  Hopkinson  bar  apparatus  used  in  this  program  is 
essentially  identical  to  that  described  by  Nicholas  (Reference  5) 
and  Bless  et  al. ,  (Reference  6).  The  essential  equations  used  to 
calculate  specimen  stress  and  strain  have  been  described  by  these 
and  other  authors.  A  brief  account  is  given  here  for  completeness. 


Strain  gauges  are  placed  equi-distant  from  the 
sample.  When  the  bar  apparatus  is  used  as  a  compression  test,  the 
initial  compressive  pulse  is  transmitted  to  the  sample.  The 
incoming  pulse  is  partially  transmitted  through  the  sample,  and 
partially  reflected.  Since  the  bars  are  of  much  larger  cross 
section  than  the  specimen,  the  boundary  condition  imposed  is  nearly 
one  of  constant  velocity,  or,  equivalently,  strain  rate. 


The  specimen  stress  is  given  by: 


o 

sp 


E 


e 


t 


(1) 


where  E  and  are  the  Young's  modulus  of  the  bar  and  the  trans¬ 

mitted  strain,  respectively,  and  A^/A^  is  the  ratio  of  bar  cross 
sectional  area  to  sample  cross  sectional  area.  The  specimen  strain 
rate  is  given  by: 
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Figure  2,  Hopkinson  Bar  Specimen  Design 
(Dimensions  in  Inches) . 
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(2) 


where  A  is  the  sample  length,  and  is  the  reflected  strain. 

These  equations  show  that  the  specimen  strain  is 
obtained  by  integrating  the  strain  pulse  reflected  from  the 
specimen,  the  specimen  stress  is  proportional  to  the  transmitted 
strain  pulse,  and  the  strain  rate  is  proportional  to  the  reflected 
strain  pulse.  The  equations  show  that  high  strairt  rates  are 
obtained  by  high  striker  velocities  or  short  sample  lengths.  The 
equations  are  not  valid  for  elastic  response  because  equilibrium  is 
not  achieved  during  the  time  needed  for  the  sample  to  "ring  up"  to 
the  bar  boundary  conditions,  as  discussed  above. 

The  above  equations  give  average,  or  engineering, 
stress  and  strain.  True  stress  and  strain  can  be  computed  from: 


a 


T 


a  ( 1+e 
sp 


(3) 


and 

e_  =  £n(l+e  )  (4) 

T  sp 

(Stress  and  strain  are  negative  in  compression.)  These  equations 
are  valid  only  when  the  stress  and  strain  are  uniform  along  the 
length  of  the  specimen. 

2. 1. 2  Test  Results 

Dynamic  stress-strain  material  characterization  tests 
were  performed  on  five  materials;  1020  steel,  C1008  steel,  HYlOO 
steel,  OFHC  copper,  and  7039-T64  aluminum.  The  strain  rate  range 
for  the  tension  tests  was  from  150  to  1800  s  The  maximum  strain 

rate  range  for  the  compression  tests  was  from  160  to  4000  s  The 

results  from  these  tests  are  summarized  in  Tables  1  and  2.  Each 
entry  in  Table  1  is  the  average  of  at  least  two  separate  tests. 

The  tables  list  the  maximum  observed  stresses.  In  cases  when 
tensile  failure  occurred  on  the  first  pulse,  these  correspond  to 
conventional  ultimate  stresses. 

The  dynamic  stress-strain  results,  from  the  Hopkinson 
bar  apparatus,  were  extremely  consistent  from  test  to  test.  Figure 
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TABLE  1 

SUMMARY  OF  TENSILE  TESTS 
a.  1020  Steel 


piax  eng  stress 
( Kbar ) * 


max  eng  strain  strain  rate 

-1 

s 


5.8 

5.9 

5.4 

4.4 


0.28 

1000 

1100 

400 

0.25 

static 

b.  QFHC  Copper 


2.6 

2.8 

2.9 

2.8 

2.5 

3.5 

1.5 


0.25 

875 

0.30 

1100 

0.30 

1100 

0.30 

1100 

0.20 

725 

0.60 

1800 

0.58 

static 

c.  Cl 008  Steel 


6.30 

6.60 

4.56 

5.00 

5.32 

4.52 

5.00 

5.05 

5.27 

3.20 

3.17 


0.30 

1000 

0.48 

1750 

0.09 

290 

0.18 

550 

0.33 

1050 

0.09 

300 

0.18 

600 

0.25 

800 

0.35 

1100 

0.23 

static 

0.23 

static 
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TABLE  1 

SUMMARY  OF  TENSILE  TESTS  (continued) 
d.  HYIOO  Steel 

max  eng  stress  max  eng  strain  strain  rate 

(Kbar)*  s"^ 


8.7 

0.05 

150 

9.9 

0.11 

350 

11.3 

0.17 

640 

11.5 

0.31 

960 

12.0 

0.30 

1100 

8.8 

0.06 

180 

10.3 

0.13 

450 

11.4 

0.23 

650 

11.7 

0.31 

930 

12.0 

0.30 

1200 

8.8 

0.23 

static 

8.8 

0.23 

static 

*  1  Kbar  =  14,504  psi 
=  100  MPA 
=  10^  dynes/cm^ 


TABLE  2 

SUMMARY  OF  COMPRESSION  TESTS 


» 


strain  rate 


max  eng  stress 


max  eng  strain 


range 


(Khar)  s 


a.  1020  Steel 


10. 

0 

0.32 

2500  - 

800 

9. 

8 

0.30 

2000  - 

500 

8. 

2 

0.23 

1200  - 

300 

8. 

5 

0.25 

1300  - 

500 

b.  OFHC  Copper 


9.0 

0.6 

3500  - 

500 

9.5 

0.6 

4000  - 

500 

7.0 

0.5 

3000  - 

500 

7.5 

0.5 

3100  - 

500 

c.  C1008  Steel 

11.0 

0.5 

1700  - 

3200 

10.0 

0.4 

800  - 

2400 

8.0 

0.25 

400  - 

2000 

d.  HYIOO  Steel 


7.5 

0.08 

100  - 

160 

9.0 

0.13 

200  - 

800 

11.0 

0.14 

200  - 

1000 

4.8 

0.08 

100  - 

500 

3.2 

0.006 

10  - 

50 

10.5 

0.16 

500  - 

12  00 

12.0 

0.18 

600  - 

1500 

e.  7039-T64  Aluminum 

8.5 

0.45 

1200  - 

3300 

7.5 

0.35 

1500  - 

2400 

6.9 

0.30 

600  - 

2000 
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3  shows  the  results  from  two  tests,  conducted  at  roughly  the  same 
strain  rate,  for  each  of  the  five  materials.  It  can  be  seen  from 
the  figure  that  the  results  are  essentially  identical  and  this  is 
typical  of  most  of  the  tests  conducted  under  this  study. 

Ouasi-static  tensile  tests  were  also  conducted  to  compare  the 
stress-strain  data  with  the  dynamic  test  data.  Table  1  summarizes 
the  quasi-static  and  tensile  SHB  data.  The  compression  SHB  data 
are  summarized  in  Table  2. 

a .  1020  steel 

Previous  SHB  data  on  1020  steel  have  been 
reported  by  Nicholas  (Reference  5)  and  Bless  et  al.,  (Reference  6). 
These  data  contain  an  unusual  amount  of  scatter,  and  they  cannot  be 
used  to  evaluate  strain  rate  sensitivity.  Additional  tests  were 
conducted  in  order  to  clarify  the  behavior  of  this  material.  The 
same  stock  was  used  as  Bless  et  al.,  (Reference  6).  Tests  were 
repeated  at  each  strain-rate  level,  spanning  400  to  1200  s  .  The 
new  data  were  obtained  based  on  the  repeatability  of  tests  at  each 
level.  The  data  for  the  dynamic  tensile  tests  are  presented  in 
terms  of  an  eyeball  curve  fit  of  the  SHB  results.  A  typical  curve 
fit  and  the  actual  data  are  shown  in  Figure  4.  The  results  from 
all  tensile  tests  conducted  on  the  1020  steel,  are  shown  in  Figure 
5.  The  flow  stress  for  a  strain  rate  of  1200  s  was  approxi¬ 
mately  15%  higher  than  the  quasi-static  flow  stress,  and  less  than 
5%  higher  than  the  flow  stress  for  a  strain  rate  of  400  s 
Hence,  there  was  no  significant  change  in  yield  stress  in  the 
dynamic  regime.  However,  the  dynamic  yield  stress  is  a  few 
percent  higher  than  the  static  yield  stress. 

Compression  tests  on  1020  steel  were  con¬ 
ducted  for  strain  rates  ranging  from  300  -  2500  s  The  strain 

rate  in  compression  tests  decreased  linearly  as  a  function  of 
time,  as  shown  in  Figure  fi.  The  results  of  the  compression  tests, 
for  three  levels  of  strain  rate,  are  presented  in  Figures  7  to  9. 
These  figures  include  the  engineering  and  true  stress  -  strain 
curves,  along  with  the  strain  rate  -  time  behavior.  They  also 
indicate  a  slight  rate  dependence. 


soo 


Figure  4.  Curve  Fit  and  Original  Data  for  a 

Dynamic  Tensile  SHE  Test  on  1020  Steel. 
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Figure  5.  Results  of  the  Dynamic  and  Static  Tensile 
Tests  on  1020  Steel. 
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Figure  6.  Strain  Rate  vs.  Time  for  a  Compression 
Test  oft  1020  Steel. 


2000 


Figure  9.  Results  Front  the  Highest  Strain  Rate  Level 
Compressive  SHB  Test  of  1020  Steel. 


b. 


QFHC  copper 

Dynamic  tensile  tests  on  the  OFHC  copper  were 
conducted  at  strain  rates  ranging  from  700  to  1800  s  The  OFHC 

copper  was  annealed,  after  machining,  at  1000°F  for  2  hours  and 
oven  cooled.  Figure  10  shows  an  example  of  the  rough  data  and 
smoothed  fit.  The  results  from  all  tensile  tests  on  OFHC  copper 
are  shown  in  Figure  11.  The  dynamic  SHB  and  static  tensile 
results  indicate  that  the  OFHC  copper  is  extremely  rate  sensitive 
and  strain  hardening. 

Compression  tests  on  OFHC  copper  were  con¬ 
ducted  for  strain  rates  ranging  from  500  to  4000  s  Stress- 

strain  and  strain  rate-time  data  for  three  levels  of  strain  rate 
are  presented  in  Figures  12  to  14.  The  observed  "overshoot"  is 
believe  to  be  an  artifact  caused  by  friction  between  the  specimen 
and  the  bars. 

c .  C1008  steel 

C1008  material  was  supplied  by  the  David 
Taylor  Naval  Ship  R&D  Center  (DTNSRDC).  Data  for  C1008  steel  were 
obtained  for  strain  rates  ranging  from  290  to  1750  s  .  Static 
tensile  tests  were  also  conducted.  The  data  are  shown  graphically 
in  Figure  15. 

The  flow  stress  for  a  strain-rate  of  1750  s  ^ 
was  approximately  three  times  the  static  flow  stress.  From  the 
SHB  tensile  test  data,  it  is  clear  that  C1008  is  an  extremely  rate 
dependent  material  even  at  low  strain-rate  levels. 

Compression  tests  on  C1008,  conducted  for 
strain  rates  ranging  from  800  to  2500  s  ,  are  presented  in 
Figures  16  to  18.  Rate  dependence  again  is  evident. 

d.  HYIOO  steel 

HYIOO  steel  alloy  was  supplied  by  the  David 
Taylor  Naval  Ship  R&D  Center.  Ten  dynamic  tensile  SHB  tests  and 
one  static  test  were  conducted.  The  results  from  different  strain 
rate  level  tensile  tests  conducted  on  the  HYIOO  steel,  are  shown 
in  Figure  19.  The  HYIOO  steel  showed  moderate  rate  dependent 
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Figure  10.  Curve  Fit  and  Original  Data  for  a 

Dynamic  Tensile  Test  on  OFHC  Copper. 
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Figure  12.  Results  From  the  Highest  Strain  Rate  Level 
Compressive  SHB  Test  on  OFHC  Copper. 
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Figure  13.  Results  From  the  Middle  Strain  Rate  Level 
Compressive  SHE  Test  on  OFHC  Copper. 
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Figure  15. 


Results  of  the  Dynamic  SHE  and  Static 
Tensile  Tests  on  C1008  Steel. 
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Figure  19 


Results  of  the  Dynamic  SHE  and  Static  Tensile 
Tests  on  HYIOO  Steel. 
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behavior.  However,  the  dynamic  flow  stress  levels  were  not 
substantially  different  than  the  quasi  static  flow  stress, 
material  also  exhibited  moderate  rate  dependency  at  low  strain 
rate  levels. 

Compression  tests  were  conducted  for  strain 
rates  ranging  from  500  to  1500  s  It  was  difficult  to  obtain 

good  data  from  the  compression  test  on  the  HYIOO  steel  because 
the  high  yield  strength  of  the  material  required  high  striker  bar 
launch  velocities.  High  striker  bar  velocity  resulted  in  noise 
in  the  bar,  which,  coupled  with  the  low  signal  (equation  2), 
made  the  strain  rate  difficult  to  interpret.  Useful  data  were 
only  obtained  in  two  of  the  eight  tests,  which  are  summarized  in 
Table  2.  The  corresponding  stress-strain  curves  and  the  strain 
rate  vs.  time  for  two  compressive  SHE  tests  are  shown  in  Figures 
20  and  21. 

e.  7039-T64  aluminum 

DTNSRDC  also  supplied  7039-T64  aluminum  alloy 
for  testing.  The  SHE  samples  were  machined  from  supplied  plate 
stock.  Table  1  lists  the  tensile  data  available,  which  span  the 
strain  rates  300  to  800  s  .  The  results  are  summarized  in 
Figure  22.  It  can  be  clearly  seen  from  the  dynamic  and  static 
tensile  tests  that  the  7039-T64  aluminum  is  a  strain  rate 
independent  material  for  strain  rates  up  to  800  s  This  result 

was  also  born  out  by  the  compression  data. 

Compression  tests  were  conducted  for  strain 
rates  ranging  from  600  to  2000  s  .  Table  2  provides  a  summary 
of  the  compression  tests.  The  results  are  presented  in  Figures 
23  to  25.  These  figures  include  the  true  stress-strain  curves 
along  with  the  strain  rate  vs.  time  behavior.  Figure  26  presents 
a  summary.  Note  that  there  is  a  tendency  for  the  true  stress  to 
decrease  slightly  for  larger  strain.  This  is  believed  to  be  a 
consequence  of  "barreling"  of  the  specimen.  The  data  beyond  the 
stress  maximum  are  probably  not  reliable. 

The  results  from  the  compression  tests  agreed 
with  the  results  from  the  dynamic  and  static  tensile  tests,  as 
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Figure  22 


Results  of  the  Dynamic  SHE  and  Static 
Tensile  Tests  on  7039-T64  Aluminum. 
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Figure  26.  The  Compressive  Stress-Strain  Behavior  of 
7039-T64  for  Three  Different  Strain  Rate 
Levels . 
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shown  in  Figure  27.  This  is  consistent  with  lack  of  rate 
dependence.  This  also  confirms  the  material  isotropy;  tests 
under  tension  and  compression  showed  the  same  behavior.  Since 
the  strain  rate  history  was  different  for  the  tensile  and 
compression  tests,  isotropy  could  not  be  directly  evaluated  for 
the  visco-plastic  materials  such  as  HYlOO  and  C1008  steel. 

2.2  NECKING  HOPKINSON  BAR  TEST 

The  results  reported  in  the  proceeding  sections  were  for 

one-dimensional  stress  state  conditions.  The  ratio  of  mean 

stress  ( o  )  to  effective  stress  (a  was  equal  to  1/3.  This 

m  eff 

stress  trajectory  limits  the  ability  to  use  Hopkinson  bar 

experiments  to  interpret  impact-induced  tensile  fracture  in  flat 

plate  impact  tests  where  the  mean  stress  to  effective  stress 

ratio  is  high  (>>  1).  Because  of  this,  available  data  generally 

show  that  the  value  of  plastic  strain  associated  with  fracture  is 

often  several  times  greater  in  a  SHB  test  than  in  a  plate  impact 

test.  Hence,  there  is  a  critical  need  to  develop  experimental 

techniques  to  obtain  dynamic  stress-strain  data  at  intermediate 

values  of  o  /a 

m  eft 

2.2.1  Analysis 

Conventional  tensile  tests,  in  which  load  and  dis¬ 
placement  are  measured,  cannot  be  used  to  determine-  constitutive 
parameters  beyond  the  point  of  maximum  load.  As  Considere 
(Reference  7)  first  observed,  an  instability  occurs  at  maximum 
lo^d.  Beyond  this  point,  the  increase  in  strength  due  to  strain 
ha  "dening  cannot  compensate  for  the  increase  in  st  "ess  due  to 
area  reduction  and  necking  occurs.  The  additional  sample 
elongation  (post  elongation)  is  greater  if  the  mat'^rial  flow 
stress  increases  with  increasing  strain  rate. 

As  the  neck  develops,  the  material  away  from  the 
neck  begins  to  unload,  and  most  of  the  plastic  deformation  takes 
place  in  a  localized  region.  To  evaluate  the  stress  state  at 
ductile  failure,  it  is  necessary  to  describe  the  deformation  in 
the  neck  region. 
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Figure  27,. 


Results  for  Static  Tensile,  Dynamic  Tensile  SHE, 
and  Dynamic  Compressive  SHE  Tests  for  7039-T64 
Aluminum. 


The  classical  work  of  Bridgman  (Reference  8)  provides 

an  approximate  method  to  obtain  the  effective  stress  and  the  strain 

at  the  neck  from  the  measurement  of  the  minimum  radius  of  the  neck 

(a)  and  the  radius  of  curvature  (R)  of  the  neck  profile.  The  axial 

strain,  e  ,  across  the  neck  is  assumed  uniform  and  equal  to: 
z 

ej  =  2  ln(d/d^)  (5) 

z  o 


where  d  is  the  current  neck  diameter  and  d  is  the  original 

o 

diameter.  The  effective  plastic  strain,  can  be  expressed  by 

the  following  relationship: 


(6) 


The  plastic  incompressibility  condition  (e^  +  =  0)  leads 

to: 


eP 

r 


(7) 


Combining  equations  (6)  and  (7),  the  effective  plastic  strain 
becomes , 


(8) 


The  effective  stress,  o  is  related  to  the  true 

ef  f 

stress,  a^,  which  is  equal  to  the  load  (as  determined  from  the  SHB 
strain  gage  data)  divided  by  the  area  of  the  minimum  cross 
section. 


a 


eff 


B  a. 


(9) 


where  the  correction  factor  B  is  given  by 
®  "  (l+2R/a)  ln(l+a/2R) 

Here  a  is  half  the  diameter  and  R  is  the  radius  of  curvature  of  a 
circle  that  osculates  the  specimen  at  the  necked  silhouette.  Since 
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B  and  are  known,  can  be  calculated  from  equation  9.  The 

effective  stress  vs.  effective  strain  curve  can  be  constructed  from 
Bridgman's  equations. 

Hancock  and  Mackenzie  (Reference  9)  conducted  tests 
on  notched  tensile  specimens  under  quasi-static  loading  conditions 
and  used  Bridgman's  analysis  to  calculate  the  stress  states  at 
ductile  failure.  The  validity  of  the  Bridgman  approximation  was 
independently  checked  by  Norris  et  al.  ,  (Reference  10).  It  was 
found  that  at  100  percent  strain,  the  error  in  flow  stress  computed 
by  the  Bridgman  analysis  was  only  10  percent. 

It  is  important  to  validate  the  extension  of  Bridgman's 
solution  for  high  rate  deformation.  For  this  purpose,  a  finite 
element  analysis  (FEA)  simulation  of  a  dynamic  tensile  test  v/as 
performed  with  the  MAGNA  (Reference  11)  program.  Attempts  were  also 
made  to  use  a  finite  difference  program,  called  STEALTH  (Reference 
12)  for  further  validation  of  Bridgman  analysis. 

a.  MAGNA 

The  first  objective  of  the  FEA  computation 
was  to  validate  the  Bridgman  analysis  for  dynamic  deformation. 

The  evaluation  of  the  neck  was  numerically  simulated.  The 
computed  contours  were  then  treated  as  input  data  for  a  Bridgman 
analysis.  The  effective  stress  that  was  inferred  from  the 
Bridgman  analysis  was  compared  with  the  actual  value  used  in  the 
code.  If  the  Bridgman  analysis  is  accurate,  the  stress  values 
will  compare  well. 

To  aid  the  calculation,  a  shallow-notched 
tensile  specimen  geometry  as  shown  in  Figure  28  was  considered. 

The  corresponding  finite  element  mesh  for  one  quarter  of  the 
specimen  is  shown  in  Figure  29.  The  number  of  4-node 
isoparametric  elements  was  440.  The  total  number  of  nodes  were 
441.  Finer  elements  were  used  around  the  minimum  cross  section 
area.  Transient  dynamic  analysis  was  carried  out  with  the  MAGNA- 
code.  Using  the  'RESTART'  facility  provided  by  MAGNA,  the 
solutions  were  carried  out  for  a  duration  of  250  microseconds 
through  several  computer  runs.  The  full  Newton-Raphson  method 
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Figure  29.  Finite  Element  Mesh  for  Shallow- 
Notched  Specimen  Geometry. 
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was  used  to  calculate  the  stiffness  matrix.  This  means  that  the 
global  stiffness  matrix  was  updated  at  each  equilibrium 
iteration.  The  solution  was  stable  and  convergence  was  achieved 
with  a  maximum  of  5-7  iterations.  The  convergence  tolerances 
were  through  displacements  and  nodal  forces.  Stringent  tolerance 
limits  were  specified  for  each  run  to  obtain  a  stable  solution 
with  the  increasing  applied  displacement.  The  full  displacement 
history  of  the  specimen  that  was  used  in  these  simulations  is 
given  in  Figure  30. 

For  Bridgman's  analysis,  contours  of  the 
notch  at  time  intervals  50,  100,  150,  200,  and  250  microseconds 
were  plotted  from  the  MAGNA  post-processor  files.  A  polynomial 
curve  was  fit  to  the  profiles  of  the  notch  at  the  five  time 
intervals.  Values  of  a  and  R  were  estimated  from  these  profiles. 
The  effective  strain  and  effective  stresses  were  calculated  using 
equations  8  and  9.  The  in  equation  9  is  readily  available 
from  the  FEA  results.  The  stress-strain  curve  obtained  through 
the  Bridgman  equations  and  the  corresponding  stress-strain  data 
used  in  the  FEA  compared  extremely  well  as  can  be  seen  from 
Figure  31. 

The  FEA  was  carried  out  for  a  maximum  strain 
of  50  percent.  Beyond  this  strain  value,  the  time  step  for  a 
stable  solution  had  to  be  very  much  smaller  and  the  time  steps 
for  stable  solution  became  unreasonably  small.  The  cost  to 
extend  this  solution  was  unacceptable  and  the  FEA  was  terminated. 
However,  the  validity  of  the  extension  of  the  Bridgman  solution 
to  high  rate  deformation  was  reasonably  justified  by  these 
limited  results. 

b.  STEALTH 

Since  the  simulation  of  a  notched  tensile 
specimen,  through  the  finite  element  code  MAGNA  turned  out  to  be 
time  consuming  and  very  expensive,  the  general  purpose,  state-of- 
the-art  finite  difference  code,  "STEALTH'  was  considered  as  an 
alternative.  The  first  objective  was  to  investigate  the  ability 
of  STEALTH  to  simulate  the  notch  tensile  test. 
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DISPLACEMENT  (mm) 


TIME  (MICROSECONDS) 


Figure  30.  Displacement  Used  in  the  Finite  Element 

Analysis  of  the  Shallow-Notched  Hopkinson 
Bar  Specimen. 


EFFECTIVE  STRESS  (KBAR> 


The  finite  element  mesh  that  was  used  in 
MAGNA  was  selected  as  a  finite  difference  grid  for  STEALTH.  A 
typical  plot  of  this  grid  from  STEALTH  is  shown  in  Figure  32. 

The  solutions  were  carried  out  for  15  microseconds.  The  material 
was  idealized  as  elastic-perf ectly  plastic.  The  pressure  (mean 
stress)  and  the  axial  stress  distributions  along  the  radial 
distance  from  the  axis  at  time  5,  10,  and  15  microseconds  of  the 
deforming  specimen  are  shown  in  Figures  33  and  34.  It  can  be 
seen  from  the  figures  that  the  distributions  which  were  smooth 
initially  become  unsmooth.  This  could  have  been  caused  by  the 
spurious  wave  reflections  from  the  free  boundaries.  The 
numerical  noise  associated  with  the  finite  difference  scheme 
could  be  another  reason.  The  pressure  and  stress  vs.  time  at  the 
axis  is  plotted  in  Figures  35  and  36.  These  plots  show  clearly 
the  spurious  numerical  noise  due  to  continuous  wave  reflections 
from  the  boundaries. 

The  undeformed  (t=0)  and  deformed  specimen 
configurations  at  t=5,  10,  and  15  are  shown  in  Figure  37.  As  the 
specimen  was  stretched  axially,  the  ratio  of  the  radius  of  a 
minimum  cross  section  to  the  radius  of  curvature  of  the  notch 
profile  increased. 

The  solutions  were  stopped  at  15  microseconds 
due  to  the  numerical  oscillations  of  the  pressure  and  stresses. 
Efforts  have  been  undertaken  to  minimize  the  effects  of  numerical 
noise  on  the  results. 

The  FEA  results  substantiated  the  extension 
of  Bridgman  analysis  to  dynamic  loading  regimes.  The  next  step 
was  to  develop  techniques  for  obtaining  the  parameters  a  and  R 
from  necking  SHB  tests.  As  mentioned  earlier,  for  this  purpose, 
a  novel  high  speed  photographic  system  was  developed  and  the 
details  are  discussed  in  the  following  section. 

2. 2. 2  Technique 

Until  now,  the  experimental  complexities  involved 
in  dynamic  tensile  tests  have  precluded  the  extension  of  the 
Bridgman  analysis  to  high  strain  rate  deformation.  To  overcome 
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Figure  32.  Grid  Used  By  the  Finite  Difference  Code  'STEALTH' 
To  Simulate  the  Shallow-Notched  SHB  Specimen. 
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Figure  33 


Mean  Stress  Distribution  Along  the  Radius 
of  the  Minimum  Cross  Section  from  STEALTH 


a.  t  =  5  ys  b.  t  =  10  ys 


c.  t  =  15  ys 


Figure  34. 


Axial  Stress  Distribution  Along  the  Radius 
of  the  Minimum  Cross  Section  from  STEALTH. 


Figure  37.  Notched  Specimen  Configuration  at  Different 
Time  Intervals  from  STEALTH. 
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this  hurdle,  a  photographic  system  using  an  LED  (light-emitting 
diode)  illuminator  and  a  35mm  rotating  drum  camera  has  been 
developed  by  Cross  et  al.,  (Reference  13). 

The  deformation  of  SHB  specimens  is  cylindrically 
symmetric,  so  back  lit  silhouette  photography  is  adequate.  Since 
back  lit  photography  is  very  efficient,  use  of  relatively  weak 
sources  is  feasible.  It  was  found  that  a  Fairchild  FLV104  LED 
(light-emitting  diode)  cooled  with  liquid  nitrogen  was 
sufficient.  The  LED  was  driven  with  pulse  currents  up  to  6amp 
with  durations  of  up  to  100ns  at  repetition  rates  of  5kHz.  The 
peak  emmision  wavelength  of  the  LED  is  670nm  which  requires  a 
film  with  an  extended  red  sensitivity.  Kodak  2479  was  selected. 

The  back-lit  photographs  were  obtained  with  a  35mm 
rotating  drum  camera.  The  camera  is  continuous  access;  framing 
was  accomplished  by  pulsing  the  illuminator.  The  drum  speed  was 
300  m/s.  The  magnification  was  about  1.2.  The  resolution  in  the 
object  plane  was  40  Ip/mm.  The  SHB  specimens  themselves  were 
about  3.2mm  in  diameter  so  that  the  specimen  diameter  could  be 
determined  to  about  one  part  in  128.  Figure  38  illustrates  a 
sample  photograph  of  a  necking  specimen. 

Synchronization  is  critical  for  successful 
correlation  of  neck  profile  measurements  with  load-displacement 
records.  Signals  from  the  two  strain  gauges  and  synch  pulses 
indicating  each  LED  flash  were  recorded  on  digital  oscilloscopes. 
In  this  way,  precise  registration  of  the  two  strain  gauge  signals 
and  the  various  photographic  images  was  possible. 

The  strain  gauge  signals  are  synchronous  with  each 
other;  however,  unlike  the  photographs,  they  are  not  "real  time" 
devices  because  of  the  finite  time  required  for  a  sound  wave  to 
propagate  from  the  sample  to  the  gauge  locations.  The  offset 
time  could  be  read  directly  from  the  strain  gauge  records,  for  it 
is  one-half  the  time  between  an  arrival  of  a  wave  traveling  to 
the  sample  and  the  arrival  of  the  corresponding  reflection. 

Prints  of  the  photographic  images  for  data  analysis 
were  conveniently  prepared  with  a  commercial  microfilm  printer. 


57 


In  order  to  calculate  R,  cross  sections  were  digitized  and  fit 
with  a  second  order  curve.  The  sample  radius,  a,  was  determined 
from  measurements  taken  from  the  photographic  prints. 

It  was  observed  that  the  uncertainty  in  the 
measurement  of  'a'  was  about  one  percent.  The  resulting  error  in 
the  strain  due  to  this  uncertainty  was  sensitive  to  the  strain 
level.  For  a  strain  of  0.01,  the  relative  error  could  be  100 
percent;  whereas,  for  a  strain  of  0.5,  the  relative  error  was 
about  4  percent.  Thus,  the  uncertainty  in  the  large  strain 
measurements  is  reasonably  small.  An  uncertainty  in  stress  of  2 
percent  results  from  an  uncertainty  in  'a'  of  1  percent;  thus, 
errors  in  'a'  do  not  cause  appreciable  errors  in  stress.  Stress 
is  also  little  affected  by  errors  in  R.  Even  for  large  defor¬ 
mations,  a  10  percent  uncertainty  in  R  causes  an  uncertainty  in 
stress  of  only  2  percent.  Thus,  significant  errors  do  not  arise 
from  uncertainties  in  the  measurements. 

2.2.3  Test  Results 

High  speed  photographs  of  necking  specimens  were 
obtained  for  each  of  the  five  metals  included  in  this  study. 
Maximum  striker  bar  velocities  were  used  in  order  to  obtain  the 
highest  possible  uniform  strain  rates  before  necking,  1000  to 
1500  s  Three  tests  were  performed  for  each  material.  The 

results  are  summarized  in  Table  3.  In  Table  3,  "uniform"  refers 
to  strains  measured  away  from  the  necking  zone.  The  term  'local' 
is  used  to  indicate  that  the  measurements  were  from  the  necking 
area.  In  some  of  these  tests,  the  onset  of  necking  was  observed 
during  the  first  tensile  pulse.  In  other  materials,  significant 
growth  of  local  necking  was  not  observed  during  the  first  pulse. 
In  Table  3,  the  strain  at  the  onset  of  necking  was  based  on 
discernable  deviations  in  the  photographs.  The  actual  onset 
could  have  occurred  slightly  before  this  time  as  can  be  seen  from 
the  strain  vs.  time  plots  reported  in  the  following  sections. 

For  C1008  steel,  OFHC  copper,  and  the  7039-T64  aluminum  the  onset 
of  necking  occurred  during  the  second  tensile  pulse.  In  these 


TABLE  3 


SUMMARY  OF  NECKING  HOPKINSON  BAR  RESULTS 


Material 

Test  # 

Uniform 
Strain  at 
Onset  of 
Necking 
(%) 

Local 

Strain 

at 

Failure 

(%) 

Loc  a  1 
Strain  at 
End  of 

1st  Pulse 
(%) 

Effective  Stress 
at  End  of 

1st  Pulse 
( Kbar ) 

1020 

HB-185 

21. 

80 

80. 

8.3 

1020 

HB-186 

24. 

96 

96. 

8.6 

1020 

HB-187 

24. 

85 

85. 

8.1 

OFHC 

HB-173 

60. 

210 

33. 

3.8 

OFHC 

HB-180 

59. 

200 

32. 

3.7 

OFHC 

HB-181 

54. 

210 

53. 

4.2 

HYIOO 

HR-241 

28. 

105 

62. 

11.7 

HYIOO 

HB-242 

36. 

104 

56. 

10.9 

HYIOO 

HB-249 

23. 

124 

66. 

12.0 

cioos 

HB-243 

52. 

134 

38. 

5.8 

C1008 

HB-244 

52. 

138 

43. 

6.0 

C1008 

HB-248 

— 

— 

-- 

— 

7039-T64 

HB-245* 

— 

— 

— 

— 

7039-T64 

HB-246 

25. 

54 

18. 

5.3 

7039-T64 

HB-247 

25. 

52 

18. 

5.4 

•k 

Poor  Photographs  made  'a'  and  'R'  calculations  impossible. 
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cases,  the  local  and  the  uniforiti  strains  are  the  same  since 
necking  has  not  been  initiated. 

1020  steel 

Three  necking  SHB  experiments  were  conducted 
on  1020  steel;  the  results  are  summarized  in  Table  3.  Local 
strain  at  the  minimum  cross  section  of  the  specimen  was 
determined  from  the  photographs.  It  was  compared  with  the 
average  strain  calculated  from  the  Hopkinson  bar  data,  as  shown 
in  Figure  39.  The  local  strain  agrees  with  the  average  strain 
until  the  onset  of  necking  at  about  70  ys.  After  the  onset  of 
necking,  the  local  strain  increases  rapidly.  The  maximum  value 
of  local  strain  is  almost  70  percent.  At  this  point,  the 
specimens  failed.  As  the  specimen  necks,  the  strain  rate 
accelerates.  For  strain  above  30  percent,  the  average  strain 
rate  is  approximately  5000  s 

Bridgman  observed  that,  for  a  number  of 
materials,  a  relationship  exists  between  the  effective  strain  and 
the  ratio  of  specimen  diameter  to  the  radius  of  curvature  of  the 
neck.  Figure  40  contains  a  plot  of  Bridgman's  results  and  the 
results  for  the  experiments  on  1020  steel.  The  dynamic  data 
falls  below  the  trend  observed  by  Bridgman  in  static  tests.  This 
is  apparently  due  to  the  material  behavior,  rather  than  inertial 
effects,  because  similar  data  for  6061-T6  aluminum  were  found  to 
lie  above  the  Bridgman  data  (Bless  et  al. ,  Reference  6). 

The  effective  stress  was  calculated  as 
described  by  equations  9  and  10.  The  effective  stress  vs. 
localized  strain  and  true  stress  vs.  localized  strain  plots,  for 
each  experiment,  are  shown  in  Figure  41.  There  is  little  change 
in  flow  stress  at  the  large  strain  and  high  strain  rates 
associated  with  the  neck  instability.  The  specimens  failed  at  an 
effective  stress  of  8.3  Kbar  and  a  localized  strain  of  85%,  as 
summarized  in  Table  3. 
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Figure  40.  Comparison  of  Bridgman's  Observed  Necking 

Behavior  and  the  Necking  Behavior  of  1020  Steel. 
(Solid  line  is  mean,  dotted  lines  are  extremes.) 
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Figure  41.  The  Effective  Stress  Calculated  from  the  Necking 
1020  Experiments. 


b. 


OFHC  copper 

Three  necking  experiments  were  conducted  on 
the  OFHC  copper,  and  the  results  are  summarized  in  Table  3. 

These  experiments  were  conducted  at  a  high  striker  bar  velocity, 
yet  necking  did  not  occur  during  the  first  pulse.  The  local 
strain  in  the  neck  is  compared  to  the  average  strain  in  Figure 
42.  The  strains  from  the  two  independent  techniques  compared 
very  well  for  each  test. 

At  the  end  of  the  first  pulse,  the  tensile 
specimens  had  reached  a  strain  of  30%  to  50%  with  no  observable 
necking.  The  difference  between  the  results  for  tests  HB-173  and 
HB-180  and  test  HB-181  is  due  to  the  higher  strain  rate  of  test 
HB-181  (1400  s  ^  as  opposed  to  800  s  ^).  During  the  second 
pulse,  necking  began  at  strains  near  60%  and  proceeded  rapidly  to 
a  localized  failure  strain  of  200%.  This  result  agrees  with 
those  of  Bauer  and  Bless  (Reference  14)  and  Fyfe  and  Rajendran 
(Reference  15)  who  also  observed  inhibition  of  necking  in  rate 
sensitive  materials.  The  stress' at  failure  could  not  be 
determined  since  SHB  strain  gauge  data  is  recorded  only  during 
the  first  tensile  pulse. 

c.  C1008  steel 

Three  necking  experiments  were  conducted  on 
the  C1008  steel,  and  the  results  are  summarized  in  Table  3. 

These  experiments  were  also  conducted  at  a  high  striker  bar 
velocity,  yet  necking  did  not  occur  during  the  first  pulse.  The 
strain  calculated  from  this  photographic  record  agreed  well  with 
the  strain  computed  from  the  SHB  gauges,  as  shown  in  Figure  43. 

At  the  end  of  the  first  pulse,  the  tensile 
specimens  had  reached  a  strain  of  40%  with  no  observable  necking. 
During  the  second  pulse,  necking  began  at  strains  greater  than 
50%  and  proceeded  rapidly  to  a  failure  strain  of  about  130  percent 
at  the  necking  section. 
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Figure  42.  Comparison  of  the  Strain  Calculated  from  the  Hopkinson  Bar 
and  from  the  Photograph  for  OFHC  Copper. 


d. 


HYIOO  steel 


Three  necking  experiments  were  conducted  on 
the  HYIOO  steel,  and  the  results  are  summarized  in  Table  3.  This 
local  strain  is  compared  with  the  average  strain  calculated  from 
the  Hopkinson  bar  data,  as  shown  in  Figure  44.  The  local  strain 
agreed  with  the  average  strain  until  the  onset  of  necking,  which 
in  this  case  occurred  at  about  25  percent  strain.  The  maximum 
value  of  local  strain  observed  is  about  65  percent. 

Figure  45  contains  a  plot  of  Bridgman's 
results  and  the  results  for  the  experiments  on  HYIOO  steel.  The 
data  from  the  HYIOO  steel  fall  within  the  scatter  of  Bridgman's 
data. 

The  effective  stress  vs.  localized  strain  and 
the  true  stress  vs.  localized  strain  plots  are  shown  in  Figure  46 
for  each  experiment.  Since  this  material  showed  significant  rate 
sensitivity,  the  flow  stress  increased  considerably  as  the  neck 
developed.  The  sharp  rise  in  the  flow  stress  is  due  to  the 
increased  strain-rates  during  necking.  The  maximum  observed 
value  of  effective  stress  was  14  Kbar  at  70  percent  localized 
strain. 

The  specimens  did  not  fail  during  the  first 
pulse.  However,  they  failed  in  the  subsequent  tensile  pulse  at 
local  strains  of  around  100  percent,  as  summarized  in  Table  3. 

The  relatively  early  onset  of  necking  in  HYIOO  steel,  compared  to 
C1008  steel,  is  consistent  with  the  observation  that  HYIOO  steel 
is  less  rate  dependent  than  C1008  steel. 

e.  7039-T64  aluminum 

Three  necking  experiments  were  conducted  on 
7039-T64  aluminum,  and  the  results  are  summarized  in  Table  3. 

The  comparison  of  local  and  average  strain  is  shown  in  Figure  47. 
The  strain  from  the  two  independent  sources  compares  very  well 
until  necking  begins. 

The  aluminum  specimens  did  not  begin  to  neck 
until  the  second  pulse  at  a  strain  of  25  percent.  Failure 


69 


a/R 


EFFECTIVE  STRAIN 


Figure  45.  Comparison  of  Bridgman's  Observed  Necking 
Behavior  (with  the  Scatter  Band)  and  the 
Necking  Behavior  of  the  HYlOO  Steel. 
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Figure  47.  Comparison  of  the  Strain  Calculated  from  the  Photographs  with 
the  Strain  Calculated  from  the  Hopkinson  Bar  for  Two  Tests 
of  7039-T64  Aluminum  at  a  Strain  rate  550  s“^. 


occurred  at  a  local  strain  of  50  percent.  Failure  appeared  to  be 
brittle.  Larger  strains  at  failure  were  observed  in  6061-T6 
aluminum  (Reference  15)  than  in  this  aluminum.  This  is  true  even 
though  the  yield  stress  of  both  materials  is  relatively 
insensitive  to  strain-rate.  One  explanation  is  that  necking  in 
'  6061-T6  is  triggered  by  void  formation  which  does  not  occur  in 

7039-T64.  Void  initiated  necking  has  been  observed  by  Bluhm 
*  (Reference  16)  in  ductile  metals  at  low  rates.  Ductile  void 

growth  is  also  known  to  occur  in  6061-T6  aluminum  (Reference  15). 

2.3  PLATE  IMPACT  EXPERIMENTS 

In  the  plate  impact  test,  a  flat  flyer  plate  is  made  to 
impact  against  a  target  plate  at  a  high  velocity.  The  flyer  and 
target  may  be  of  the  same  or  different  material.  Compressive 
stresses  are  produced  and  transmitted  immediately  from  the  plane 
of  impact  to  the  adjacent  stress  free  areas  of  the  material  in 
the  form  of  a  stress  pulse.  Many  discussions  of  planar  impact 
loading  are  available  (References  17  and  18). 

Plate  impact  test  provide  a  loading  path  that  is  very 
different  from  conventional  SHB  compression  or  tension  tests. 

The  deformation  is  that  of  one-dimensional  strain,  and  the  mean 
stress  is  generally  very  high  compared  to  that  is  SHB  tests. 
Strain  rates  are  10^  s  ^  or  higher.  The  material  undergoes 
compression  immediately  followed  by  tension.  Thus,  plate  impact 
experiments  are  essential  for  calibrating  and  validating  high 
strain  rate  material  models  that  aspire  to  general  applicability. 
Specifically,  plate  impact  data  may  be  interpreted  to  infer 
■«  compression  and  tensile  yield  strengths  and  failure  parameters. 

2.3.1  Techniques 

T 

The  experimental  techniques  for  plate  impact 
experiments  at  the  University  of  Dayton  are  described  by  Bless 
(Reference  19)  and  Bless  et  al.  ,  (Reference  6).  Readers  who  are 
unfamiliar  with  this  technique  should  look  at  Figure  48  which 
provides  an  (x,t)  diagram.  The  origin  represents  the  instant  of 
impact  on  the  target  surface.  Shock  waves  propagate  into  the 
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target  and  into  the  flyer  plate.  Shock  waves  reflect  from  free 
surfaces  as  release  (tensile)  waves.  In  the  example  diagrammed, 
these  lead  to  a  spall  at  point  SP.  Shock  waves  emanate  from  the 
spall  plane  as  the  stress  relaxes  to  zero.  The  velocity  of  the 
rear  of  the  target  is  measured  with  a  VISAR  (Reference  20).  In 
some  experiments,  a  transparent  window  is  placed  behind  the 
target  in  order  to  suppress  release  waves  and  prevent  spall. 

2.3.2  Yield  Strength  and  Spall  Threshold 

The  plate  impact  tests  were  conducted  with  three 
objectives:  (1)  determination  of  Hugoniot  Elastic  Limit,  (2) 

determination  of  the  unloading  path  from  the  free  surface 
velocity  history,  and  (3)  determination  of  the  threshold 
conditions  associated  with  onset  of  spall  fracture. 

Impact  induces  an  elastic  shock  and  a  plastic  shock 
in  the  target.  The  amplitude  of  the  elastic  shock  is  The 

Hugoniot  elastic  limit,  ,  is  the  maximum  stress  for  one- 

dimensional  elastic  wave  propagation.  This  stress  level  is  a 
material  property,  and  above  this  level  the  material  flows 
plastically.  The  stress  can  be  determined  from  the 

experimentally  obtained  free  surface  velocity  of  the  target  that 
corresponds  to  the  elastic  shock, 
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where  is  the  elastic  sound  speed.  The  high  strain  rate  yield 
strength  can  then  be  calculated  from  the  relationship: 
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where  G  is  the  shear  modulus  and  K  is  the  bulk  modulus. 

In  conventional  elastic  plastic  theory,  the  release 
process  is  initially  elastic.  Release  waves  travel  at  the 
elastic  wavespeed.  The  elastic  deviatoric  stress  possible  in 
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release  is  2Y.  However,  not  all  materials  display  this  behavior. 
Brittle  materials  in  particular  are  sometimes  observed  to  undergo 
loss  of  shear  strength  behind  the  initial  shock.  In  order  to 
determine  the  release  path  from  data,  the  material  is  modeled  by 
a  numerical  method.  The  constitutive  descriptions  are  varied 
until  the  observed  free  surface  velocity  decay  associated  with 
release  wave  arrivals  is  reproduced. 

The  "spall  stress"  a^,  is  defined  as  the  highest 
tensile  stress  experienced  by  the  target  prior  to  spall.  The 
spall  stress  is  often  computed  from 

a  =  l/2p  C,  AV  (13) 

S  Lj  s 

(see  Figure  49  for  AV  definition).  However,  wave  propagation 

s 

calculations  show  that  this  is  incorrect  if  spall  occurs  abruptly 
because  shock  waves  emanating  from  the  spall  plane  will  overtake 
release  waves. 

Equation  (13)  also  is  based  on  the  assumption  that 
all  of  the  release  characteristics  which  arrive  at  the  rear 
surface  before  the  spall  signal  propagate  at  speed  (Reference 
19).  In  fact,  only  the  lead  characteristic  has  this  speed.  This 
assumption  causes  equation  (13)  to  overestimate  by  possibly  as 
much  as  ten  percent. 

2.3.3  Test  Results 

Plate  impact  experiments  were  conducted  on  C1008 
and  HYIOO  steels,  and  7039-T64  aluminum,  see  Table  4.  A  few 
tests  were  conducted  on  1020  steel  to  supplement  previous  data 
(References  6,19)  and  support  the  numerical  analysis. 

Preliminary  experiments  were  also  conducted  with  a  double  flyer 
plate  technique  to  observe  void  formation  and  recompaction. 

Target-samples  (77mm  diameter  and  6mm  thickness) 
were  fabricated  for  each  of  the  metals  to  be  tested.  The  target 
plates  were  mounted  in  the  fixture  which  was  attached  to  the 
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TABLE  4 


PLATE  IMPACT  TESTS  ON  C1008,  HYlOO,  and  7039-T64 


Shot 

No. 

Velocity 

(m/sec) 

Target* 

Projectile 

Visible 

Spall 

C1008  Steel 

695 

157 

6  mm 

3  mm 

C1008 

1020 

696 

373 

6  mm 

3  mm 

C1008 

1020 

678 

178 

6  mm 

3  nun 

part. 

C1008 

1020 

685 

134 

6  mm 

3  mm 

no 

C1008 

1020 

687 

134 

6  mm 

3  mm 

no 

C1008 

1020 

679 

192 

6  mm 

3  mm  1020 

n.a* 

C1008W 

3  mm  Lexan 

677 

200 

6  mm 

3  mm  1020 

n.a. 

C1008W 

1.6  mm  WC 

HYlOO  Steel 

694 

203 

6  mm 

3  mm 

part. 

HYlOO 

1020 

680 

355 

6  mm 

3  mm 

yes 

HYlOO 

1020 

681 

180 

6  mm 

3  mm 

no 

HYlOO 

1020 

686 

200 

6  mm 

3  mm 

part. 

HYlOO 

1020 

688 

196 

6  mm 

3  mm  1020 

n.  a. 

HYIOOW 

3  mm  Lexan 

689 

187 

6  mm 

3  mm  1020 

n.  a. 

HYIOOW 

1.6  mm  WC 

690 

186 

6  mm 

3  mm  1020 

n.a. 

HYIOOW 

3  mm  Lexan 

7039-T64 

Aluminum 

693 

180 

6  mm 

3  mm 

yes 

7039-T64 

1100 

682 

375 

6  mm 

3  mm 

yes 

7039-T64 

1020  St 

683 

187 

6  mm 

3  mm 

part. 

7039-T64 

1100  A1 

684 

162 

6  mm 

3  mm 

no 

7039-T64 

1100  A1 

691 

200 

6  mm 

3  mm  1100 

n.a. 

7039-T64  W 

3  mm  Lexan 

692 

200 

6  mm 

3  mm 

n.a. 

7039-T64  W 

1.66  mm  WC 

*  W  denotes  Si02  window  behind  the  target 
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muzzle  of  the  50mm  compressed  air  gun  as  explained  earlier. 
Projectile  flyer  plates  were  mostly  3mm  thick. 

In  a  few  tests,  the  flyer  plates  were  mounted  on 
tungsten  carbide  (WC)  backing  plates  so  that  the  initial  shock 
loading  would  be  followed  by  reshock  instead  of  release.  Using 
both  release  and  recompression  data  permits  a  check  for  relax¬ 
ation  of  shear  stress  behind  the  shock  waves.  Many  targets  were 
backed  with  fused  silica  windows  6mm  thick  and  25mm  in  diameter. 
The  windows  prevent  release  waves  from  originating  at  the  target 
rear  surface,  and  thus  suppress  spall.  This  set-up  allows 
observation  of  the  complete  release  process,  which  might 
otherwise  be  blocked  by  spall  fracture.  The  windows  remained 
transparent  at  least  until  the  initial  shock  wave  had  passed  all 
the  way  through  them,  which  gave  a  recording  time  of  about  one 
microsecond.  The  windows  were  depolished,  silvered,  and  bonded 
to  the  targets  with  thin  film  epoxy.  A  stop  was  placed  on  the 
VISAR  operator  to  block  the  spectral  reflections  from  the  window 
rear  surface.  The  VISAR  data  collected  were  used  to  assist  in 
determining  the  elastic  yield,  spall  threshold,  and  spall 
strength. 

a.  1020  steel 

Experiments  on  1020  steel  have  been  chiefly 
motivated  by  a  need  to  fully  describe  spall  fracture.  There  are 
two  approaches  for  deriving  spall  criteria  in  plate  impact  tests, 
and  both  have  been  utilized.  The  first  is  to  vary  the  flyer 
plate  thickness  and  velocity  and  determine  the  threshold  impact 
conditions  for  target  damage.  It  is  usually  found  that  thinner 
flyer  plates  require  higher  impact  velocities  to  produce  damage. 
The  threshold  criteria  can  be  interpreted  as  the  combinations  of 
stress  level  and  duration  just  sufficient  to  cause  fracture. 

The  second  approach  is  to  measure  transmitted 
wave  profiles  after  damage  has  occurred.  Target  sectioning  and 
wave  profile  data  determine  the  time  and  location  where  spall 
first  occurred. 
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The  two  approaches  to  spall  characterization 
are  sampling  two  rather  different  types  of  behavior.  Threshold 
studies  generally  employ  relatively  long  duration  tensile  pulses. 
Transmitted  wave  techniques  operate  on  a  much  shorter  time  scale, 
and  spall  separation  may  be  limited  by  time  dependent  processes. 
Experiments  were  performed  with  1020  steel  of  hardness  62  to  75 
on  the  Rockwell  B  scale.  In  most  experiments,  a  conventional 
flyer  plate  was  used.  However,  in  some  tests,  the  flyer  plate 
had  a  step  along  a  diameter  on  its  rear  surface.  With  this 
technique  each  impact  was  equivalent  to  two  experiments  with 
conventional  flyer  plates.  In  other  shots,  tapered  flyer  plates 
were  used.  The  motive  here  was  not  so  much  to  obtain  damage 
gradients,  as  in  Seaman  et  al.  ,  (Reference  21),  but  rather  to 
efficiently  define  the  spall  threshold  in  regions  in  which  the 
most  important  parameter  is  flyer  plate  thickness.  This  is  the 
case  for  thin  flyer  plates;  for  thick  plates,  the  damage  is 
principally  velocity-dependent.  Use  of  wedge  flyers  has  enabled 
us  to  obtain  data  for  shorter  loading  pulses  than  are  normally 
reported  in  spall  experiments.  The  wedge  technique  is  described 
in  more  detail  by  Bless  et  al. ,  (Reference  6). 

The  microstructure  of  the  steel  consisted  of 
equiaxial  grains  whose  typical  diameter  was  30  ym.  The  flyer 
plates  were  backed  with  PMMA  or  air,  except  the  wedges,  which 
were  backed  with  mating  teflon  wedges.  Free  surface  velocity  was 
measured  in  the  center  of  the  target  disc.  Complete  spall  was 
defined  as  a  void  content  of  over  50  percent  on  the  spall  plane. 
Partial  spall  was  defined  as  visible  (at  400X)  damage  less  than 
complete  spall.  Results  reported  below  are  mainly  based  on 
additional  analyses  of  data  in  Reference  6. 

The  o  data  for  1020  steel  were 

nhii-t 

inconclusive,  in  spite  of  a  large  number  of  shots.  Figure  50 
shows  a  typical  velocity  profile.  The  HEL  is  not  well  defined, 
even  though  the  target  was  9mm  thick  and  the  shock  rise  was 
resolved.  SHB  results  indicate  a  yield  strength  of  about  4.5  Kbar 
and  work  hardening  to  7  Kbar;  7  Kbar  corresponds  to  Kbar 

and  u„„.  =  0.05  mm/ys.  This  arrival  is  not  evident  in  Figure  50; 

rltj  J_» 
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Figure  50.  VISAR  Data  for  1020  Steel,  showing  poorly 
Defined  HEL.  (Target  was  9mni  Thick,  Flyer 
Plate  was  4inm.) 
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nor  is  it  apparent  on  other  VISAR  traces.  However,  reference  22 
reports  an  HEL  of  11  Kbar  for  12.  7min  thick  targets  of  1020  steel 
very  similar  to  our  material. 

The  scatter  in  HEL  data  is  probably  due  to  the 
"precursor  decay  phenomena"  that  has  been  observed  in  iron 
(Reference  23).  Data  for  Armco  iron  from  several  sources  (e.g.  , 
References  23  and  24)  indicate  that  in  5mm  thick  targets,  the 
amplitude  of  the  elastic  precursor  is  about  50  percent  higher  than 
its  equilibrium  value. 

It  appears  that  the  "spall  signal"  in  1020 
steel  is  generated  by  void  nucleation.  The  evidence  is  provided 
by  shot  529  which  was  at  400°C.  The  VISAR  data  is  shown  in  Figure 
51.  A  spall  reverberation  is  clear.  After  sectioning  and 
polishing,  this  target  appeared  undamaged.  An  electron  microscope 
V7as  necessary  to  reveal  micron  sized  pores  as  shown  in  Figure  52. 
These  were  widely  distributed  about  the  incipient  spall  plane. 

The  critical  velocity  for  spall  development 
can  be  read  from  Figure  53.  The  flyer  plates  used  in  these 
experiments  were  backed  with  PMMA.  The  threshold  for  partial 
spall  was  about  200  m/s  for  2mm  flyer  plates,  and  between  251  and 
279m/s  for  1mm  thick  flyer  plates.  Values  of  spall  stress  from 
equation  (13)  were  only  available  from  three  shots  at  impact 
velocities  about  500m/s.  The  range  was  26  to  33  Kbar.  The  most 
negative  stress  corresponding  to  incipient  spall  with  2mm  thick 
flyers  is  about  22  Kbar.  Reference  22  found  that  for  3.2mm  flyer 
plates,  incipient  spall  occurred  at  15  to  17  Kbar.  In  Figure  53, 
this  corresponds  to  an  impact  velocity  of  about  150m/s. 

b.  OFHC  copper 

No  additional  plate  impact  experiments  were 
performed  with  copper  beyond  those  reported  by  Bless  et  al. , 
(Reference  6).  However,  additional  analysis  of  their  data  was 
carried  out.  Figure  54  summarizes  the  observation  of  spall 
failure  in  this  material.  The  spall  stress,  ,  computed  from 
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51. 


VISAR  Data  for  Shot  529,  5iTiin  1020  Target 
2min  Flyer.  (Target  Temperature  was  400 °C 


52.  Electron  Microscope  Photograph  of 
Micron  Sized  Pores  in  1020  Target, 
Shot  529,  lOOOx. 
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Figure  53.  Critical  Velocity  for  Spall  Development 
of  1020  Steel.  (Solid  circles  denote 
complete  system;  half-filled  circles, 
partial  spall;  open  circles,  no  spall. 
Flag  shows  widely  separated  micron-size 
voids . ) 
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equation  (13)  was  between  12  and  13  Kbar  for  all  of  the  copper 
records. 

c.  C1008  steel 

The  evaluation  of  threshold  conditions 
required  several  shots  with  different  velocities.  Five  tests  were 
conducted  to  determine  the  spall  threshold  velocity.  For  a  better 
description  of  the  shocked  material's  strength,  two  more  tests 
were  performed  with  the  flyer  plates  mounted  on  WC/Lexan  backing 
plates  as  explained  in  the  preceding  paragraph. 

The  VISAR  was  used  on  all  shots  except  695, 
the  VISAR  results  are  shown  in  Figure  55.  A  complete  release  wave 
profile  was  obtained  in  shot  679  because  spall  was  suppressed  by 
the  silica  window  and  Lexan  backing  plate  behind  the  flyer  plate. 
The  release  profile  clearly  shows  the  arrival  of  two  elastic  waves 
in  agreement  with  conventional  elastic-plastic  theory. 

Most  of  the  VISAR  traces  show  a  ramp  shock 
wave.  Other  impact  tests  on  mild  steel  and  iron  have  not  shown 
this  anomaly,  although  the  data  here  are  at  slightly  lower  impact 
velocities  than  those  reported  elsewhere.  This  feature  of  the 
waves  is  so  unusual  that  it  should  be  checked  in  further 
experiments. 

Four  good  observations  of  the  elastic 
precursor  were  obtained.  The  value  of  the  HEL  varied  from  11  to 
12  Kbar.  The  average  value  was  11.3  Kbar.  However,  as  pointed 
out  in  the  section  on  1020  steel,  the  true  value  of  the  HEL  may  be 
50  percent  less  than  this. 

The  spall  threshold  was  bounded  by  shots  at 

157m/s  and  178m/s.  At  the  lower  velocity,  the  voids  were 

connected  by  a  continuous  crack.  The  threshold  velocity  for  3mm 

thick  flyer  plates  is  thus  168  ±  lOm/s.  The  corresponding  spall 

stress  is  31  ±  2  Kbar.  This  stress  should  be  used  to  predict 

spall  when  damage  accumulation  is  not  taken  into  account  and  the 

tensile  pulse  duration  is  about  lys.  The  value  of  a  calculated 

s 

from  equation  (13)  was  18  Kbar  for  shots  678,  685  and  687. 
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The  C1008  spall  threshold  is  consistent  with 
the  1020  spall  data.  This  is  shown  in  Figure  56.  The  curve  in 
Figure  56  extended  to  3min  thick  flyer  plates  also  shows  the  SHE 
data  for  1020  steel.  The  point  is  plotted  at  the  impact  velocity 
that  would  have  produced  the  same  mean  stress  as  occurred  in  the 
SHE  specimen  at  failure 

cj^  =  1/3  (1  +  3  Jin  (1  +-^))  <7eff 

The  agreement  between  these  data  sets  shows  that  the  Tuler-Eutcher 
criterion  used  with  a  mean  stress  formulation  is  capable  of 
describing  failure  in  both  of  these  mild  steels  over  stress  states 
from  one-dimensional  stress  to  one-dimensional  strain.  The  mean 
stress  corresponding  to  cr^  in  the  Tuler-Butcher  equation  is 

a  =  cr  -  |y  (15) 

m  o  3 

d.  HYIOO  steel 

Five  tests  were  carried  out  on  HYIOO  steel 
target  plates.  VISAR  data  were  obtained  in  all  the  HYIOO  shots, 
as  shown  in  Figure  57.  In  shot  690,  the  data  were  difficult  to 
interpret,  and  since  that  shot  was  nearly  the  same  as  688,  the 
data  from  the  latter  were  used  for  analysis. 

The  shock  waves  in  the  silica-backed  targets 
do  not  appear  steady.  We  believe  that  this  is  probably  an 
artifact  of  the  distortion  introduced  by  the  shock  wave  traversing 
the  window.  The  release  or  reshock  in  these  experiments  appears 
sharp  and  unambiguous. 

The  HEL  was  well  resolved  in  shots  680  and 
681.  In  both,  the  value  was  15.5  ±  .3  Kbar.  The  release  profiles 
in  shots  680  and  681  clearly  indicate  elastic  arrivals  of 
amplitude  similar  to  the  initial  HEL,  in  agreement  with  elastic 
plastic  theory.  However,  the  reshock  in  shot  689  is  preceded  by  a 
release  wave  that  travels  at  an  elastic  velocity,  indicating  that 
an  air  gap  was  present  between  the  two  flyer  plates. 
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Figure  57.  VISAR  Data  from  Impacts  on  HYIOO  Steel  Plates 


Unfortunately,  the  presence  of  the  air  gap  eliminates  the 
possibility  of  measuring  the  speed  of  the  reshock. 

The  spall  threshold  for  HYIOO  is  bracketed  by 

shot  681,  which  resulted  in  isolated  voids,  and  shot  686,  which 

caused  a  spall  crack.  The  threshold  velocity  is  thus  190  ±  10 

mm/ys,  corresponding  to  a  simple  spall  stress  of  35  ±  2  Kbar.  The 

value  of  a  from  equation  (18)  seemed  to  vary  with  impact 
s 

velocity:  in  shot  681,  a  =27  Kbar;  in  shot  686,  a  =28  Kbar, 

s  s 

and  in  shot  680,  a  =31  Kbar. 

s 

e.  7039-T64  aluminum 

Three  shots  were  made  with  various  velocities 
to  determine  the  spall  threshold  of  7039-T64  aluminum.  Two 
additional  shots  were  done  to  obtain  release  and  reshock  wave 
profiles.  VISAR  data  were  obtained  as  shown  in  Figure  58. 

Three  values  of  the  HEL  were  recorded;  8  Kbar 
in  shot  683,  8.4  Kbar  in  shot  691,  and  6.8  Kbar  in  shot  692.  The 
weighted  mean  value  is  7.7  Kbar.  As  was  the  case  for  steel,  the 
data  from  shot  683  show  separated  elastic  and  plastic  wave 
arrivals,  and  the  first  elastic  wave  has  about  the  same  amplitude 
as  the  HEL.  Note  that  the  release  in  shot  683  is  essentially 
complete,  since  no  spall  occurred.  Spall  also  did  not  occur  in 
shot  684.  Again,  the  release  is  nearly  completed  and  a  ringing 
period  corresponding  to  a  double  wave  transit  through  the  target 
is  present. 

Comparing  shots  691  and  692  (Lexan  and  WC 
hacked  flyer  plates)  is  difficult  because  the  first  arrival  in 
shot  692  is  not  clear.  The  most  reliable  time  comparison  is 
between  the  midpoint  of  the  plastic  shock  and  the  midpoint  of  the 
arrival  from  the  flyer/backing  plate  interface.  This  is  0.94  ys 
in  shot  691  (release)  and  0.87  ys  in  shot  692  (compression). 

Given  the  width  of  the  wavefront,  the  time  difference  of  70  ns 
between  these  two  transit  times  is  probably  not  significant.  The 
expected  delay  time  between  the  arrival  of  the  plastic  wave  and 
the  arrival  of  the  wave  from  the  rear  surface  of  the  aluminum 
flyer  plate  is  0.94  ys,  ignoring  the  relatively  weak  elastic  wave 
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Figure  58. 


VISAR  Data  from  Five  Impacts  on  7039-T64 
Aluminum  Targets. 


in  the  flyer  plate.  The  expected  difference  if  the  release  is 
wholly  plastic  in  the  target  is  1.13  ys.  Clearly  both  release  and 
recompression  waves  are  traveling  with  elastic  wavespeeds  in  the 
target,  indicating  either  substantial  relaxation  of  shear  stress 
or  hardening  of  the  material  behind  the  shock. 

The  spall  threshold  of  this  material  was  not 
well  determined.  Shot  684  suffered  no  damage.  However,  shot  682, 
which  was  meant  to  provide  an  upper  bound,  was  inadvertantly  shot 
with  a  steel  flyer  plate,  which  produced  a  much  higher  stress  than 
intended  (and  caused  complete  spall).  Thus,  only  a  lower  bound 
for  the  simple  spall  stress  is  presently  available:  13.5  Kbar. 

The  value  of  (equation  13)  was  13  Kbar  in  shot  682  (complete 
spall)  and  17  Kbar  in  shot  683  (partial  spall). 

2.3.4  Experiments  with  BeO 

The  properties  of  porous  polycrystalline  BeO  under 
shock  compression  below  and  just  above  the  HEL  were  determined. 

The  release  and  spall  behavior  were  inferred  from  analysis  of  the 
VISAR  data.  Data  for  very  high  pressures  and  material  identical 
to  that  used  here  are  available  in  Marsh  (Reference  25). 

a.  Material 

The  BeO  was  purchased  from  Ceradyne  Inc.  as 
Ceralloy  418.  It  was  supplied  as  discs  40mm  in  diameter  and 
thickness  4,  8,  and  12mm.  Eighteen  specimens  were  obtained.  They 

3 

were  ordered  with  a  minimum  density  of  2.985  ±  .005  g/cm  .  The 

3 

crystal  density  of  BeO  is  3.008  g/cm  . 

Sample  density,  P  was  measured  by  the 
Archimedes  method  of  liquid  immersion.  During  immersion,  the 
specimens  were  observed  to  gain  weight,  indicating  water 
penetration.  After  the  initial  measurements,  the  specimens  were 
boiled  in  deionized  water  for  5  minutes,  cooled,  and  left  immersed 
for  24  hours.  Immersion  weights  were  then  remeasured. 


If  we  denote  the  density  calculated  from  the 

initial  immersed  weight  as  ,  and  the  second  as  and  the  volumes 

of  solid  ceramic,  open  pores,  and  closed  pores  as  V  ,  V  ,  and  V  , 

s  o  c 

the  apparent  porosity  is  defined  as 


n 


V 

c 


(Vg  +  V^) 


(16) 


The  total  porosity  is 


(V^  +  V^) 
c  o 

(V  +  V  +  V  ) 
s  c  o 


(17) 


Values  for  n  are  listed  in  Table  5. 

For  thick  disc  samples,  n  equaled  n^.  For 
the  8mm  samples,  n^  was  typically  20%  greater  than  n.  For  4mm 
samples,  n^  was  about  three  times  n.  The  values  of  n^  were 
similar  for  all  samples  -  about  2%.  However,  the  data  were  not 
consistent  with  either  of  these  two  assumptions:  (1)  there  is  a 
constant  porosity  associated  with  the  surfaces,  or  (2)  all 
samples  have  the  same  values  of  n^,  and  the  liquid  diffused 
similar  distances  into  each  specimen,  resulting  in  differing 
values  of  n. 

We  have  resolved  this  indeterminacy  by 
assuming  that  in  each  specimen,  the  value  of  n  represents  the 
interior  porosity.  When  n^  exceeds  n,  this  is  assumed  to  be 
caused  by  surface  porosity.  For  shock  wave  analysis,  therefore, 
n  was  considered  to  be  the  relevant  porosity. 

Sound  speed  was  measured  by  a  pulse-echo 
technique.  The  results  are  also  listed  in  Table  5.  It  was  found 
that  sound  speed  correlated  with  n,  as  shown  in  Figure  59.  The 
correlation  substantiates  the  assumption  that  n  is  a  measure  of 
the  interior  porosity  of  the  specimens.  It  is  shown  in  Figure  59 
that  the  extrapolated  zero-porosity  sound  speed  is  11.93  mm/ys. 
This  is  in  good  agreement  with  literature  values:  11.99  (Marsh, 
Reference  25),  12.05  (Soga,  Reference  26),  and  11.54  (Gust  and 
Royce,  Reference  27).  The  dashed  line  in  the  figure  is  an 
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TABLE  5 

POROSITY  OF  BeO  SPECIMENS 


Specimen 

number 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


Thickness 

(mm) 

12 

12 

8 

8 

8 

8 

8 

8 

8 

8 

4 

4 

4 

4 

4 


Porosity 
n  (%) 

2.7 

1.9 

1.4 

1.6 

1.4 

1.4 

1.4 

1.2 

1.5 

1.1 

0.6 

0.6 

0.6 

0.6 

0.6 


Sound  Speed 
CL  (mm/ys) 

11.94 

11.72 

11.85 

11.43 

11.85 

11.47 

11.64 

11.85 

11.85 

11.64 

11.74 

11.74 

11.74 

11.74 

11.74 
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Figure  59.  Observed  Variation  of  BeO  Sound 
Speed  with  Porosity. 
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extrapolation  to  data  of  Marsh,  which  were  for  much  greater 
porosity. 

b.  Shock  Compression  Experiments 

Table  6  lists  the  experiments  conducted  with 
BeO  targets.  Flyer  plates  were  annealed  OFHC  copper  or  annealed 
1100  aluminum. 

Preliminary  tests  with  aluminum  oxide  ceramic 
targets  (Reference  28)  demonstrated  that  VISAR  data  could  not  be 
obtained  from  a  ceramic  free  surface.  The  intensity  and  contrast 
appeared  to  drastically  degrade  when  the  plastic  shock  arrived  at 
the  rear  target  surface.  It  was  not  even  possible  to  measure  the 
magnitude  of  the  peak  particle  velocity.  Therefore,  several 
targets  were  backed  with  PMMA  windows  for  use  with  the  VISAR. 

Thin  film  epoxy  was  used  to  bond  the  windows.  The  side  adjacent 
to  the  target  was  depolished  and  vacuum-coated  with  aluminum. 

The  rear  surface  was  cut  to  a  7“  angle  to  that  the  specular 
reflection  which  it  generated  did  not  enter  the  VISAR.  Window 
corrections  published  by  Barker  (Reference  20)  were  used  in  the 
VISAR  data  reduction.  Other  targets  were  backed  by  a  copper 
disc.  VISAR  measurements  were  made  on  the  rear  surface  of  the 
copper. 

Figure  60  illustrates  an  x,t  diagram  and  o,u 
diagram  for  impacts  with  PMMA  backing  plates.  The  example  shown 
is  for  shots  728  and  729,  in  which  the  flyer  plate  was  aluminum. 
Point  1  denotes  the  impact.  This  point  also  denotes  the  material 
behind  the  shock  waves  emanating  from  the  point  referred  to  as 
state  1.  The  pressure  in  state  1  is  below  the  HEL  of  the 
ceramic,  so  there  is  only  a  single  elastic  shock  wave  generated 
in  the  BeO.  The  HEL  of  the  aluminum  is  so  low  that  the  elastic 
wave  in  that  material  is  ignored?  only  the  plastic  wave  is  shown. 
Point  2  represents  the  arrival  of  the  ceramic  shock  at  the 
window;  this  is  the  arrival  detected  by  the  VISAR.  The  pressure 
at  this  point  is  reduced  almost  to  zero,  and  the  velocity  is 
nearly  twice  the  velocity  in  state  1.  The  shock  in  the  aluminum 
reflects  from  the  copper  at  point  3.  It  is  at  a  higher  stress 
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Solution  Diagram  for  Shots  with  Aluminum 
Flyer  Plate  and  PMMA  Windows. 


than  state  1.  The  reflected  shock  enters  the  BeO  at  point  4. 

The  stress  is  increased  once  again,  as  shown  in  the  o,u  diagram. 
State  5  results  from  the  intersection  of  the  shock  waves  from 
points  2  and  4.  Note  that  the  state  5  stress  is  also  greater 
than  in  state  1.  Hence,  this  target  does  not  spall.  The  time  of 
the  arrival  at  point  6  has  no  direct  physical  interpretation, 
since  it  depends  on  the  speeds  of  the  release  wave  from  4  to  5 
and  the  shock  wave  from  2  to  5. 

Figure  61  shows  x,t  and  o,u  diagrams  for 
shots  731  and  732  for  copper  flyer  plates  above  the  HEL  of  the 
BeO.  The  elastic  wave  in  the  BeO  is  shown  as  1'.  Again,  the 
elastic  wave  in  the  flyer  plate  is  so  weak  that  it  can  be 
ignored.  In  this  case,  state  3  is  at  a  lower  stress  than  state 
1.  State  4  is  a  release  from  state  1.  Finally,  state  5  becomes 
tensile,  and  spall  may  be  expected  to  occur  at  point  5.  The 
VISAR  observes  points  2',  2,  and  the  spall  signal  at  point  6. 

The  HEL  (point  1')  is  calculated  by  assuming  a  symmetric  BeO 
Hugoniot.  The  plastic  stress,  point  2,  is  calculated  by  using 
the  observed  arrival  time  of  point  2  to  calculate  the  shock 
velocity  which  corresponds  to  the  slope  of  the  segment  from 
points  1'  to  1  on  the  ct,u  plot,  and  by  assuming  that  the  release 
path  is  elastic.  The  arrival  time  observed  for  state  2  indicated 
that  the  release  path  from  state  1  to  state  2  is  elastic.  This 
is  consistent  with  Asay,  Chhabildas,  Dandekar  (Reference  29)  and 
Grady  (Reference  30),  who  reported  that  even  in  elastic-isotropic 
materials,  the  initial  release  is  elastic. 

Figure  62  illustrates  the  solution  diagrams 
for  shots  with  a  copper  plate  behind  the  target.  These  are 
relatively  simple.  Observation  of  the  free  surface  arrival  at 
points  3  and  5  can  be  used  with  the  known  Hugoniot  of  copper  to 
calculate  points  2  and  4  on  the  BeO  Hugoniot.  The  slope  of  the 
measured  Hugoniot  gives  the  shock  impedance  of  the  BeO: 
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Solution  Diagram  for  Shots  with  Copper 
Flyer  Plate  and  PMMA  Windows. 


4 


Cu 


Figure 


62.  The  Solution  Diagrams  for  Shots  with  Copper 

Flyer  Plate  on  Target  with  Copper  Back  Plate. 
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Figure  63  illustrates  sample  VISAR  data  from  a  copper  back  plate. 
The  two  shock  states  are  well  defined. 

Figure  64  presents  the  o,e  Hugoniot 
calculated  from  the  experiments.  The  impedance  in  the  elastic 
region  is  351  Kbar/mm/us. 

Yaziv  (Reference  28  and  31)  showed  that  the 
shock  properties  of  fully  dense  BeO  can  be  inferred  from  these 
data  by  using  the  Steinberg  model  (Reference  32)  for  the  effect 
of  porosity  on  the  Hel.  According  to  this  model,  when  a  porous 
material  of  initial  specific  volume  V'^  is  shocked  to  its  HEL 


V  -  V'  =  V  - 

O  HEL  OO  HEL 


a  V 
HEL  o 


(19) 


where  V  refers  to  the  fully  dense  solid. 

OO 


It  follows  that 


V  =  V 
HEL  OO 


(^'hel> 
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(20) 


Steinberg  based  his  model  on  aluminum  oxide 

data  from  Gust  and  Royce  (Reference  27).  In  order  to  assess 

its  generality,  we  also  applied  the  model  to  porous  tungsten 

from  data  Linde  (Reference  33),  Dandekar  (Reference  39),  and 

McQueen  (Reference  18).  Using  data  from  these  three  sources, 

AV  =  V  -  was  computed  as  shown  in  the  table  below. 

OO  HEL 


TABLE  7 

EFFECT  OF  POROSITY  ON  HEL 


P 

o 

a 

HEL 
( Kbar ) 

^L 

AV 

(cm^/g ) 

Linde 

14.  3 

9.6 

3.  67 

. 00035 

Dandekar 

15.  29 

14.3 

4.26 

.00039 

McQueen 

19.2 

37.  5 

5.  22 

. 00037 
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The  average  value  was  .00035.  The  small  deviation  about  the 
average  demonstrates  the  accuracy  of  the  Steinberg  model.  For 
BeO,  V  -  V„__  =  0.  0074  cmVg. 

O  nhjJj 

Figure  65  compares  the  results  of  solid  HELs 
computed  by  this  technique  with  literature  data.  The  solid  HEL 
is  calculated  to  be  96  Kbar. 

The  point  labelled  S  is  the  value  inferred  by 

Steinberg  from  the  data  of  Gust.  Steinberg  obtained  this  value 

by  averaging  data  from  Gust  and  Royce  which  was  obtained  for  four 

different  target  thicknesses  -  3,  6.4,  8.9,  and  12.8mm.  We 

believe  that  the  variation  of  o  ,  with  target  thickness  is 

caused  by  precursor  decay.  If  only  the  value  from  the  3mm  thick 

target  is  used,  the  predicted  solid  o  is  within  the 

HEL 

experimental  uncertainty  of  our  value. 


c.  Spall  Measurements  in  BeO 


The  relationship  between  the  free  surface 
velocity  and  the  spall  strength  is  affected  by  the  presence  of  a 
window  (Yaziv,  Reference  28  and  31).  Figure  66  illustrates  the 
general  forms  of  the  window/target  interface  velocity.  In  the 
case  of  no  spall,  the  trace  will  have  form  (a).  In  the  case  of 
spall,  the  form  will  be  as  in  (b).  The  spall  strength  will  be 


0 
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2^^BeO 


+  7  \ 

mMA' 


Au 


(21) 


(where  and  are,  respectively,  the  shock  impedances  of 

the  ceramic  and  the  window)  if  shock  waves  generated  at  the  spall 
surface  do  not  overtake  release  waves  originating  from  the  flyer 
rear  surface.  When  the  shock  compression  remains  elastic,  this 
should  be  the  case,  and  spall  strength  can  be  directly  computed 
from  the  spall  signal  by  the  above  equation. 


Figure  67  shows  an  x, t  diagram  for  spall  in 
an  elastically  compressed  solid.  A  release  fan  is  generated  in 
the  metal  flyer  plate.  However,  in  the  ceramic,  all  of  the 
release  characteristics  propagate  at  the  same  Lagrangian 
velocity.  The  release  characteristic  arrives  at  the  observed 
interface  at  point  R.  The  wave  that  reflects  back  into  the 
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Figure  67.  Lagrangian  x,t  Diagram  for  Spall  in 
Ceramic  Target. 
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target  at  point  R  is  a  shock.  Spall  first  occurs  at  the  point  S. 

A  shock  is  generated  that  arrives  at  the  window  interface  at  time 

t  .  The  location  of  point  S  can  be  computed  from  t  ,  since  the 
s  s 

wave  velocities  are  known  to  all  be  equal  to  C  .  The  wave  that 

reflects  from  the  window  at  time  t  is  a  release  wave. 

s 

Spall  was  investigated  in  two  experiments 
that  provided  especially  clear  records:  shots  731  and  732. 

Figure  68  shows  the  interface  velocity  history  from  shot  731, 
where  there  is  a  weak  plastic  wave  also  present.  The  spall 
strength  in  732  is  1.8  Kbar.  The  spall  strength  in  shot  731, 
which  was  compressed  above  the  HEL,  is  negligible.  These  data 
show  that  BeO  loses  its  spall  strength  above  the  HEL. 

Yaziv  (Reference  28  and  31)  showed  that  a 
regional  spall  develops  behind  the  spall  signal  that  originates  at 
point  S  in  Figure  67.  The  spall  is  quenched  by  the  arrival  of  the 
shock  wave  that  originates  at  point  R.  The  trapped  stress  waves 
"ring",  as  shown  in  Figure  66.  The  distance  can  be  computed 
from  the  ringing  period  of  the  spall  signal.  At  ,  while  x  can  be 

5  S 

computed  from  t  .  In  this  way,  the  VISAR  record  can  be  used  to 
s 

compute  the  distance  Ax  =  x^  -  x^  in  Figure  67.  In  shot  731, 

Ax  =  2mm,  In  shot  732,  Ax  =  1.75mm.  Both  of  these  values  are  in 
quantitative  agreement  with  the  spall  region  hypothesis. 

Although  the  spall  data  are  limited,  they 
point  to  a  consistent  picture  of  spall  behavior.  The  spall 
strength  for  elastically  stressed  material  is  about  1.8  Kbar. 

This  value  diminishes  above  the  HEL.  Lack  of  spall  strength 
above  the  HEL  has  also  been  observed  by  Munson  and  Lawrence 
(Reference  35)  in  alumina.  In  the  elastic  region,  spall  occurs 
not  in  a  plane,  but  over  a  region  that  can  be  predicted  by  simple 
characteristics  analysis. 

2.4  NEW  EXPERIMENTAL  TECHNIQUES 

In  this  section  we  present  developments  of  improved 
techniques  for  high  strain  rate  characterization.  However,  these 
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techniques  have  not  yet  been  used  in  the  material  modeling 
efforts. 

2.4.1  Double  Flyer  Plate  Technique 

Preliminary  experiments  were  conducted  to  show 
feasibility  of  a  powerful  new  technique  to  study  the  properties 
of  shock  damaged  materials.  Progress  on  this  task  was  also 
reported  by  Yaziv  and  Bless  (Reference  36). 

Two  flyer  plates  are  separated  by  a  small  gap 
(Figure  69).  The  first  plate  causes  spall  in  the  target,  as  in  a 
conventional  experiment.  The  second  plate  is  of  higher  shock 
impedance.  Impact  of  the  second  plate  closes  the  fracture.  The 
advantage  of  this  technique  is  that  it  allows  determination  of 
the  properties  of  the  damaged  material  in  the  spall  region. 

Preliminary  experiments  have  been  conducted  to 
demonstrate  feasibility.  A  2-inch  (50mm)  diameter  compressed  gas 
gun  was  used  to  launch  50mm  diameter  flyer  plates.  Target  rear 
surface  motion  was  measured  with  the  VISAR.  The  wave  inter¬ 
actions  in  the  target  and  flyer  plates  were  anlaysed  with  the 
SWAP  method  of  characteristics  code.  The  target  was  fully 
annealed  OFHC  copper.  The  first  flyer  plate  was  2024-T3 
aluminum,  and  the  second  was  copper  or  1020  steel. 

Five  shots  were  completed  that  demonstrated 
fracture/recompaction  in  copper.  The  velocity  range  was  250  to 
400  m/s.  Targets  were  0.16,  4mm  thick  and  flyer  plates  were  2mm 
thick.  Shots  719  and  720  were  typical;  the  velocities  were  365 
and  388  m/s,  respectively.  Figure  70  shows  the  x,t  diagram 
generated  by  the  SWAP  code  for  shot  719,  In  this  shot  and  in 
shot  720,  only  half  of  the  aluminum  first  flyer  was  backed  by  a 
copper  second  flyer;  the  other  half  was  air  backed.  Figure  71 
shows  the  appearance  of  the  targets  after  impact.  The  two  halves 
of  the  target  display  different  damage,  as  predicted  by  the  SWAP 
simulation.  The  single  flyer  plate  impact  produced  partial  spall 
(a  region  containing  porosity),  while  no  damage  is  readily 
visible  opposite  the  two-flyer  plate  impact. 
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Figure  69.  Experiment  Configuration  for  Double 
Flyer  Impact  Experiment. 
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Figure  71. 


target  target 
FLYERS  #720  #719 


Shot  719  and  720  After  Impact. 
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The  free  surface  velocity  record  from  shot  719  is 
shown  in  Figure  72.  In  addition  to  the  conventional  "spall 
signal"  (point  3)  there  is  an  arrival  from  the  second  flyer 
impact  (point  4)  that  has  "leaked"  through  the  damaged  region, 
and  a  shock  arrival  (point  6)  caused  by  final  closure  of  the 
spall  plane. 

2. 4. 2  Cranz-Schardin  Camera 

A  Cranz-Schardin  camera  was  developed  in  order  to 
obtain  high  resolution  sequential  photographs  of  projectiles. 

The  intended  application  was  photography  of  Taylor-impact 
specimens. 

The  basic  configuration  of  the  Cranz-Schardin 
camera  is  illustrated  in  Figure  73.  An  array  of  short-duration 
light  sources,  generally  closely  spaced,  produces  light  which 
fills  a  main  lens.  The  main  lens  is  designed  to  focus  the  source 
array  onto  a  lens  array.  A  small  diameter  camera  lens  is 
positioned  at  each  of  the  source  images  so  that  as  each  source  is 
pulsed,  a  single  camera  lens  receives  the  light  that  it  produces. 
The  light  produced  by  all  the  sources  overlaps  in  the  region 
close  to  the  main  lens.  The  subject  is  located  close  to  the  main 
lens,  and  the  small  camera  (objective)  lenses  are  designed  to 
form  images  of  the  object  on  a  photographic  film.  As  each  light 
source  pulses,  a  distinct  image  of  the  object  is  recorded.  By 
firing  the  light  sources  in  sequence,  the  time  evolution  of  an 
object  may  be  studied  via  back-lit  photography. 

Conventional  Cranz-Schardin  cameras  employ  an  array 
of  spark  gaps  as  sources.  Instead  of  using  spark  sources,  the 
University's  camera  incorporates  commercially-available  Light- 
Emitting  Diodes  (RED-1000).  These  sources  are  very  small, 
require  relatively  little  power,  are  semiconductor-compatible, 
and  should  require  no  maintenance  for  thousands  of  photographs. 

a.  Light-Emitting  Diodes  (LED) 

A  variety  of  commercially-available  LED's 
have  been  investigated  as  possible  sources  for  the  Cranz-Schardin 
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Figure  72.  Free  Surface  Velocity  for  Shot  719. 
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Figure  73. 


Cranz-Schardin  Camera  Operation.  Two  of  array  of 
twenty  sources  and  objective  lenses  are  shown. 
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camera.  It  was  found  that  the  manufacturers'  specifications  for 
LED  time-response  were  unreliable.  Several  manufacturers'  data 
sheets  list  a  variety  of  response  times  for  LED's,  and  we  found 
that  invariably  these  were  very  conservative,  i.e. ,  the  quoted 
time  response  was  very  much  slower  than  that  which  was  actually 
observed  when  the  LED's  were  operated  with  a  current-driving 
circuit  which  was  capable  of  producing  very  short  pulses.  This 
is  what  one  would  expect  of  a  solid  state  source  which  produces 
light  by  electron-hole  recombination.  Outside  of  a  small 
resistance  (a  few  ohms)  and  a  small  capacitance  (a  few  hundred  pF 
at  most)  there  is  no  basic  physical  process  which  should 
practically  limit  the  turn-on  and  turn-off  times  of  solid  state 
emitters.  We  chose  to  test  only  those  commercially  available 
LED's  which  had  the  highest  brightness  in  the  forward  direction. 
Among  those  tested,  no  correlation  was  found  between  the 
specified  time  response  and  that  of  the  measured  response  when  an 
LED  was  used  with  the  proper  circuitry. 

Some  light  emitting  diodes  exhibit  a 
saturation  in  the  light  output  intensity  as  a  function  of 
current,  as  can  be  seen  in  manufacturers'  specification  sheets. 
LED's  exhibiting  this  behavior  were  avoided. 

In  pulsed  operation,  the  magnitude  of  the 
current  pulse  is  at  least  a  factor  of  100  higher  than  the  DC 
current  values  at  which  the  performance  of  an  LED  is  usually 
specified.  Therefore,  pulsed  performance  tests  of  candidate 
LED's  were  necessary  in  order  to  identify  those  which  emitted  the 
highest  energy  light  flash.  It  was  found  that  the  specified  DC 
light  output  values  were  often  very  misleading;  in  the  pulsed 
mode,  "high  intensity"  LED's  may  perform  much  poorer  than  "lower 
intensity"  ones. 

Infrared  laser  diodes  were  also  considered  as 
sources.  In  general,  the  total  amount  of  light  which  can  be 
produced  by  an  Infrared  Laser  diode  is  an  order  of  magnitude 
greater  than  available  from  LED's.  However,  extended  red 
sensitivity  films  have  such  superb  speed  and  resolution  (compared 


to  infrared  films)  that  the  practical  advantages  of  laser  diodes 
are  greatly  diminished.  In  the  extreme  case  of  a  very  rapidly 
moving  subject,  however,  a  nanosecond  duration  infrared  laser 
diode  would  be  the  only  usable  source. 

Two  types  of  emitter-film  combinations  were 
tried  in  the  camera;  (1)  a  long-wavelength  LED  whose  peak 
emission  wavelength  is  about  660  nm  and  Kodak  2485  film,  and 
(2)  a  shorter  wavelength  LED  whose  emission  wavelength  is  about 
630  nm  and  Polaroid  film.  The  total  light  energy  output  from  the 
latter  LED  was  very  much  smaller  than  from  the  former  but  the 
much  higher  sensitivity  of  Polaroid  film  made  up  for  the  decrease 
in  light  emission.  The  use  of  Polaroid  film  is  very  convenient 
but  a  penalty  is  paid  in  terms  of  resolving  power,  since  Polaroid 
film  can  resolve  only  about  25-30  Ip/mm  whereas  the  Kodak  film 
can  resolve  about  50. 

b.  Camera  Description 

The  final  Cranz-Schard in  camera  design  used 
an  array  of  FED-1000  LED's  arranged  in  a  4  x  5  pattern  on  0.283 
inch  centers.  The  LED's  are  individually  driven  by  a  set  of 
GA301  (4)  silicon  controlled  rectifiers  fired  in  sequence.  The 
sequencer  consists  of  a  set  of  three  type  74164  IC  shift 
registers  driven  by  an  external  pulse  generator.  Initially,  a 
logical  "one"  is  fed  into  the  first  bit  of  the  first  74164  IC. 

All  the  other  bits  of  the  other  chips  are  logical  "zeroes".  As 
pulses  are  fed  into  the  first  shift  register,  the  "one"  shifts 
over  from  bit  to  bit  and  as  each  transition  is  made,  a  trigger 
pulse  is  generated  which  is  fed  into  one  of  the  GA301  SCR's.  The 
shift  registers  have  been  tested  up  to  a  shifting  rate  of  23  MHz. 
Usually,  however,  the  use  of  such  a  high  framing  rate  is 
unnecessary.  In  the  experiments  described  below,  no  rate  higher 
than  500  kHz  has  been  used. 

The  discharge  of  each  SCR  dumps  a  0.01-mFd 
capacitor  charged  to  between  40  and  60  volts.  The  maximum 
holdoff  of  the  GA301  SCR  is  80  volts  but  it  has  been  found  that 
little  increase  in  the  total  light  output  from  an  LED  occurs 


above  a  60  volt  capacitor  charging  voltage.  The  light  output 
pulse  from  each  LED  has  a  FWHM  of  60  ns  indicating  an  effective 
resistance  for  the  LED  of  approximately  3  ohms  at  a  peak  current 
of  nearly  10  amp.  A  smaller  image  size,  say  1.25-cm  diameter, 
would  have  allowed  the  capacitor  to  be  reduced  in  size  by  about  a 
factor  of  four  and  yielded  a  light  pulse  with  a  considerably 
shorter  pulse  duration.  However,  in  the  Taylor  impact  test,  the 
subject  velocity  was  generally  less  than  500  m/s.  Maximum 
subject  spatial  resolution  was  of  primary  importance,  so  temporal 
resolution  was  sacrificed.  For  subject  velocity  of  500  m/s,  the 
motion  blur  due  to  the  exposure  time  was  insignificant. 

The  main  lens  of  the  camera  is  a  simple 
biconvex  glass  lens  of  100mm  diameter  and  167mm  focal  length. 

The  separation  between  the  main  lens  and  the  plane  of  the  LED's 
is  31.8  cm.  The  LED's  are  therefore  imaged  at  a  magnification  of 
about  3  onto  a  set  of  small  camera  lenses  spaced  on  20mm  centers. 
The  camera  lenses  are  19mm  diameter  achromats  of  190mm  focal 
length.  The  separation  between  the  main  lens  and  the  plane  of 
the  camera  lenses  is  approximately  102  cm.  The  camera  lenses 
were  adjusted  to  have  a  field  of  view  approximately  7.0  cm  in 
diameter  located  about  19  cm  from  the  main  lens.  The  shadow 
image  of  any  test  object  is  therefore  recorded  at  a  magnification 
of  about  0.36.  The  camera  is  relatively  insensitive  to  alignment 
and  precision  optical  adjustments  are  not  necessary.  A  calibra¬ 
tion  target  was  always  located  in  the  field  of  view  so  that  each 
individual  image  could  be  calibrated  for  absolute  magnification 
and  parallax. 

Kodak  2485  4  inch  by  5  inch  sheet  film  was 
used  as  the  recording  medium  and  it  was  developed  under  the 
standard  conditions  specified  by  the  manufacturer;  D-19 
developer,  95°F  for  2.5  minutes.  The  first  was  a  prototype  for 
proof  of  principle  testing.  It  was  used  for  the  Taylor  impact 
study  described  in  the  following  section.  It  is  pictured  in 
Figure  74a.  The  second  camera  was  a  self-contained  unit  that  was 


delivered  to  AFATL/DLJW.  It  was  housed  in  a  modified  4x5 
camera  body.  Figure  74b  shows  this  unit. 

In  the  delivered  camera,  the  spacing  between 
the  LED's  was  0.250  inches  to  allow  for  a  longer  distance  between 
the  main  lens  and  the  subject  volume.  The  main  lenses  diameter 
and  focal  length  were  unchanged.  The  camera  lenses  were  changed 
to  16mm  diameter  achromats  of  124mm  focal  length  to  allow  the 
working  distance  of  the  camera  to  be  increased. 

c.  Photographic  Results 

The  camera  has  been  employed  in  two  types  of 
impact  studies,  both  of  which  have  involved  small  steel  rod 
impactors.  The  first  of  these,  described  as  a  Taylor  Impact 
Test,  involves  launching  a  small  metal  rod  onto  the  surface  of  a 
hardened  steel  impact  plate  and  studying  the  deformation  of  the 
rod  as  time  proceeds.  In  our  experiments,  the  rod  was  of 
annealed  Armco  iron  (nearly  impurity  free  iron)  launched  at 
velocities  between  135  and  360  m/s  from  a  0.3-inch  smooth  bore 
barrel  either  by  means  of  compressed  helium  gas  or  by  gunpowder 
(above  180  m/s).  The  rods  were  0.3  inch  in  diameter  and  one 
inch  long.  Independent  of  the  projectile  velocity,  the  framing 
rate  employed  was  500  kHz. 

A  triggering  system  located  close  to  the 
projectile  impact  point  was  required  in  order  to  run  the  camera 
relatively  fast.  The  triggering  pulse  to  initiate  the  firing  of 
the  sequence  of  LED's  was  generated  by  the  interruption  of  a  beam 
from  a  Fairchild  Type  FPE104  IR  LED.  This  LED  produces  a  very 
high  brightness,  low-divergence  beam  which  was  detected  by  a 
Motorola  type  MDR500  photodiode.  The  FPE104  and  the  MDR500  were 
positioned  so  that  the  IR  light  beam  was  interrupted  by  the 
passage  of  the  projectile  approximately  5mm  in  front  of  the 
impact  plate.  Figure  75  illustrates  a  sequence  of  several  images 
obtained  at  a  projectile  velocity  of  344  m/s.  (All  of  the  frames 
recorded  by  the  camera  are  not  illustrated. ) 

A  saw  blade  with  20  teeth  per  inch  was 
positioned  in  the  same  vertical  plane  as  the  projectile.  It 


appears  at  the  bottom  of  every  image.  The  distance  between  the 
various  teeth  of  the  blade  were  measured  with  a  microscope  and 
compared  to  those  measured  from  a  final  print  of  an  enlargement 
of  the  4x5  Cranz-Schardin  negative.  The  lateral  magnifications 
of  the  camera  depends  upon  the  cosines  of  the  angles  that  the  LED 
optic  axes  make  with  the  optic  axis  of  the  main  lens.  The 
maximum  value  of  the  angle  is  only  0.0491  radians  and  the 
correction  factor  would  be  only  0.999.  Measurements  made  on 
printed  images  verified  that  the  f rame-to-f rame  differences  in 
magnification  were  too  small  to  resolve.  For  most  data  analyses, 
the  2.5  cm  diameter  images  were  magnified  by  a  factor  of  about  8 
during  the  printing  process.  By  measuring  all  the  points  of  the 
saw  blade  and  comparing  the  digitized  values  with  the  known 
values,  an  average  value  for  the  Cranz-Schardin  magnification 
could  be  determined.  A  typical  magnification  was  0.3770.  The 
standard  deviation  of  the  magnifications  determined  for  the 
various  saw  tooth  intervals  was  0.0001. 

The  2485  film  employed  to  record  these  images 
is  reported  to  be  capable  of  resolving  50  Ip/mm.  Since  the 
object  image  ratio  of  the  camera  was  3:1,  the  resolution  in  the 
object  plane  was  about  16.7  Ip/mm  or  0.006  cm  (0.0024  inch).  At 
a  projectile  velocity  of  450  m/s  and  an  exposure  time  of  60  ns, 
the  projectiles  were  capable  of  moving  only  0.0027  cm  during  the 
recording.  This  amount  of  motion  blur  is  completely  negligible 
because  it  cannot  be  resolved  on  the  film. 

The  second  type  of  experiment  conducted  using 
this  camera  was  for  the  purpose  of  studying  crack  formation  and 
propagation  in  glass.  Steel  rods  were  fired  into  a  reservoir 
containing  fine  aluminum  oxide  powder.  The  reservoir  external 
dimensions  were  2.5  x  2.5  x  1.25  cm.  It  was  constructed  of  1.6mm 
thick  PMMA.  The  open  top  of  the  reservoir  was  capped  by  a  piece 
of  plate  glass  of  the  same  size  as  the  reservoir  and  the  two  were 
clamped  together  to  ensure  a  firm  contact  between  the  powder  and 
the  glass.  This  assembly  was  positioned  inside  the  Cranz- 
Schardin  camera  so  that  the  view  was  of  the  1x1  inch  face  of 
the  glass.  Steel  rods  0.3  inch  x  1. 0  inch  were  fired  into  the 


powder  reservoir  just  below  the  interface  between  the  glass  and 
the  powder.  When  a  rod  impacted  the  powder,  the  powder  was 
driven  ahead  of  the  rod  and  the  high  pressures  generated  by  this 
process  eventually  cracked  the  glass.  The  objective  of  the 
experiment  was  to  determine  exactly  where  and  when  the  glass 
cracked  (relative  to  the  position  off  the  rod  tip).  After  each 
experiment,  the  glass  was  always  fragmented  into  many  small 
pieces.  However,  the  photographs  illustrate  that  not  many  cracks 
are  actually  formed  during  the  impact  itself. 

Figure  76  illustrates  four  frames  recorded 
during  an  experiment.  The  projectile  velocity  was  431  m/s  and 
the  camera  framing  rate  was  225  kHz.  Small  cracks  were  always 
observed  to  form  in  the  portion  of  the  glass  nearest  the  impact 
point.  These  small  cracks  are  on  the  extreme  right  of  the  image 
of  Figure  76.  As  the  projectile  velocity  was  increased,  first 
one  and  then  two  cracks  such  as  are  illustrated  in  the  figure 
were  formed.  The  first  crack  always  propagated  upwards  through 
the  glass  and  terminated  at  the  base  of  a  screw  which  was  used  to 
hold  the  glass  and  reservoir  together.  The  propagation  velocity 
of  the  crack  through  the  glass  was  approximately  1690  m/s  as 
measured  during  these  experiments.  Figure  76  also  illustrates 
the  images  of  the  leading  edge  of  an  air  shock  which  was 
generated  by  the  impact.  This  curved  edge  is  on  the  far  left  in 
each  of  the  images  and  can  be  seen  to  propagate  further  to  the 
left  as  time  proceeds.  We  have  computed  the  velocity  of 
propagation  of  a  shock  in  air. 

2.4.3  Taylor  Impact  Tests 

Taylor  impact  tests  provide  a  technique  to  measure 
compressive  flow  stress  at  high  strain  rates  (References  37,  38). 
The  objective  of  the  University's  test  program  was  to  show  that 
the  Cranz-Schardin  camera  could  be  used  for  obtaining  profile 
data  during  the  impact  process.  A  secondary  objective  was  to 
obtain  deformation  contours  for  analysis  by  personnel  at 
AFATL/DLJW. 


Fiaure  76.  Four  Frames  from  Test  Shot  into 
Ceramic/Glass  Target. 
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The  tests  were  conducted  on  7.62min  diameter  2.54mm 
steel  projectiles  launched  from  a  1.5mm  long  7.62mm  diameter 
smooth  bore  barrel.  The  propellant  was  either  compressed  helium 
gas  or  gunpowder.  The  gas  was  used  for  velocities  below  180  m/s, 
the  gunpowder  at  velocities  above  this.  The  projectile  velocity 
was  determined  by  measuring  the  time  interval  between  the  inter¬ 
ruption  of  two  He-Ne  beams  which  illuminated  MRD500  PIN  photo¬ 
diodes. 

The  projectiles  were  formed  from  three  types  of 
steel:  (1)  common  mild  steel  (type  unspecified),  (2)  Armco  iron 

which  was  annealed  in  air,  and  (3)  vacuum-annealed  Armco  iron. 

The  Armco  iron  specimens  were  annealed  at  1700°F 
for  one  hour,  and  cooled  at  a  rate  of  250°F  per  hour.  The  first 
two  projectile  types  were  used  to  generate  a  loading  curve  for 
the  gun  and  to  conduct  tests  for  the  Cranz-Schard in  camera. 

The  steel  or  Armco  iron  rods  impacted  normally  upon 
a  massive  impact  block  which  was  machined  from  BD30  steel  and 
hardened  to  Rockwell  C53-57.  Figure  77  shows  a  photograph  of  one 
of  these  blocks.  At  first,  the  blocks  were  made  flat  across  the 
whole  front  surface.  It  was  found  that  the  images  produced  by 
the  Cranz-Schardin  camera  were  affected  by  the  extreme  edges  of 
the  block.  Experience  showed  that  aiming  accuracy  was  extremely 
good,  so  that  in  later  tests  the  block  was  bevelled  away  from  the 
impact  site. 

The  impact  yaw  could  be  detected  in  the  Cranz- 
Schardin  camera  images  and  from  the  condition  of  the  projectile 
recovered  after  the  experiment  was  over.  On  occasion,  particu¬ 
larly  at  low  velocities,  the  projectile  elastically  rebounded 
from  the  impact  plate,  out  of  the  camera,  and  onto  the  laboratory 
floor.  More  frequently,  it  was  found  to  lie  inside  the  Cranz- 
Schardin  camera  housing. 

The  deformation  of  the  rod  during  impact  was 
detected  via  the  Cranz-Schardin  camera  images.  Data  analysis 
began  with  an  examination  to  determine  which  of  the  20  frames 
should  oe  enlarged  for  further  quantitative  study.  The  selected 
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Figure  77. 
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Photograph  of  Taylor  Impact  Target. 
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images  from  each  experiment  were  enlarged  to  4  x  5  inch  prints 
and  anlayzed  using  a  Talos  digitizer.  A  computer  program  was 
developed  specifically  for  this  task.  The  methodology  of  the 
image  analysis  included  several  steps:  (1)  A  small  saw  blade  was 
located  in  the  field  of  view  of  the  camera  and  was  recorded  on 
each  image.  The  positions  of  the  teeth  tips  were  digitized  and 
compared  to  the  known  distances  between  the  teeth  as  measured 
with  a  traveling  microscope.  This  comparison  allowed  the  calcu- 
tion  of  the  average  and  variance  of  the  magnification  of  each 
image.  (2)  Two  points  defining  the  direction  of  the  surface  of 
the  impact  plate  and  two  points  defining  the  rear  end  of  the  rod 
were  digitized.  This  allowed  the  calculation  of  the  length  of 
the  rod  in  the  particular  image  and  defined  the  origin  of  the 
coordinate  system  of  the  image.  The  origin  was  defined  as  the 
intersection  of  the  rod  axis  and  the  shadow  of  the  impact  plate. 
(3)  A  grid  was  laid  down  on  the  image  and  pairs  of  points 
diametrically  opposite  one  another  were  digitized.  The 
computation  of  the  distance  between  these  points  allowed  the 
calculation  of  the  instantaneous  rod  diameter  as  a  function  of 
length  along  the  rod. 

A  total  of  38  projectile  firings  were  performed. 
Fifteen  were  for  the  purpose  of  generating  a  loading  curve  for 
the  gun  for  both  helium  and  gunpowder  propellants.  An  additional 
11  shots  were  made  while  the  operating  procedures  for  the  Cranz- 
Schardin  camera  were  being  developed.  The  final  shots  were  with 
projectiles  made  of  vacuum-annealed  Armco  iron  furnished  by 
AFATL. 

The  framing  rate  necessary  for  the  Cranz-Schardin 
camera  was  estimated.  A  typical  projectile  velocity  was  taken  as 
240  m/s.  The  projectile  was  assumed  to  collapse  to  75  percent  of 
its  original  length  during  the  impact  and  to  decelerate 
uniformly.  Thus  it  stops  in  0.6  cm  with  a  mean  velocity  during 
the  slow  down  of  126  m/s.  The  time  over  which  the  projectile 
stops  is  50.6  ys.  Therefore,  the  camera  can  record  all  of  the 
impact  in  20  frames  if  the  interframe  time  is  2.5  ms.  This 


corresponds  to  a  framing  rate  of  400  kHz.  The  rate  used  in  all 
data  tests  was  500  kHz  to  account  for  nonuniform  deceleration. 

Difficulty  was  initially  encountered  in  synchro¬ 
nizing  the  Cranz-Schardin  camera  properly  with  the  impact.  The 
last  photodiode  trigger  signal  from  the  velocity  station  was  used 
to  trigger  the  camera.  A  signal  delay  generator  was  inserted  in 
the  line  between  the  last  photodiode  and  the  trigger  input  to  the 
camera.  From  the  anticipated  projectile  velocity  and  the  known 
distance  between  the  last  timing  station  and  the  impact  plate,  it 
was  possible  to  calculate  the  time  delay  and  to  trigger  the 
camera  at  an  appropriate  time.  However,  this  turned  out  to  be  an 
unreliable  technique.  It  was  possible  to  set  the  time  interval 
of  the  delay  generator  accurately  enough  to  run  the  camera  at  a 
suitable  framing  rate.  The  chief  problem  was  the  inherent 
variations  in  the  velocity  of  the  gun. 

Precise  trigger  delays  are  usually  achieved  by 
installing  an  up-down  counter  on  the  range.  While  this  approach 
would  have  been  acceptable  for  the  experiments  conducted  at  the 
university,  where  such  a  counter  was  available,  this  solution  was 
not  acceptable  for  AFATL,  who  may  not  have  this  specialized 
device  available.  Therefore,  it  was  decided  that  another  type  of 
triggering  system  should  be  developed. 

Larger  and  larger  variations  in  the  projectile 
velocity  can  be  tolerated  by  the  camera  system  as  the  triggering 
point  is  moved  closer  to  the  impact  plate.  Accordingly,  a  third 
He-Ne  laser-Pin  diode  station  was  installed  on  the  range.  The 
beam  of  this  He-Ne  laser  was  aligned  so  that  it  crossed  the  tra¬ 
jectory  of  the  projectile  about  2.5  cm  in  front  of  the  impact 
plate.  In  this  configuration,  very  large  projectile  velocity 
variations  were  tolerable  because  the  total  time  from  when  the 
laser  beam  was  interrupted  to  when  the  projectile  struck  the 
impact  plate  was  so  short.  However,  this  technique  proved  to  be 
unacceptable  because  the  scattering  of  the  laser  light  from  the 
dust  in  the  air  caused  an  intense  streak  to  be  recorded  on  the 
film,  partially  overwriting  some  of  the  Cranz-Schardin  images. 


Nevertheless,  several  test  films  were  made  showing  that  trig¬ 
gering  the  camera  so  close  to  the  impact  point  worked  very  well, 
except  for  the  streak  on  the  film.  The  solution  that  provides  an 
adequate  trigger  was  use  of  an  infrared  diode  and  detector,  as 
described  in  the  previous  section. 

The  most  complete  set  of  experiments  were  those 
conducted  with  annealed  Armco  iron.  The  test  matrix  is  given  in 
Table  8.  There  exists  a  critical  velocity  above  which  radial 
cracks  form  along  the  perimeter  of  the  impact  face  of  the  rod. 
Taylor  impact  tests  are  normally  conducted  below  this  critical 
velocity.  The  critical  velocity  was  determined  experimentally 
from  shots  0094  to  0102.  In  shot  94,  the  velocity  was  344  m/s, 
and  severe  radial  cracking  occurred.  The  launch  velocity  was 
lowered  step-by-step  (shots  0095  to  0099).  At  225  m/s,  cracks 
formed,  but  at  198  m/s  (shots  0099  to  0100)  no  radial  cracks  were 
formed.  In  order  to  determine  the  threshold  velocity  to  greater 
accuracy,  the  launch  velocity  was  slowly  raised  until  a  region 
where  crack  formation  was  problematical  was  found.  This  region 
was  demonstrated  to  lie  around  216  m/s,  since  in  shot  0101  there 
was  no  crack  formation  at  a  velocity  of  217  m/s  but  in  shot  102  a 
crack  formed  at  215  m/s.  It  was  observed  that  a  small  chip  of 
iron  was  always  released  when  a  crack  formed.  If  the  front  edge 
of  the  sample  was  not  chipped,  there  never  was  a  crack  formed. 
End-on  and  side-on  views  of  the  impacted  rods  for  the  shots  0094 
to  0105  are  illustrated  in  Figure  78.  Figure  79  show  a  dynamic 
profile  obtained  with  the  Cranz-Schard i n  camera. 

As  illustrated  in  the  figures,  at  high  velocities  a 
debris  cloud  was  produced.  The  debris  cloud  was  unexpected  and 
tended  to  obscure  the  end  of  the  rod  near  the  plate.  Three 
additional  firings  were  conducted  at  lower  velocities  to  try  to 
obtain  images  in  which  there  was  no  cloud.  Success  was  obtained 
for  shots  0103  to  0105.  Thus,  the  velocity  at  which  the  cloud 
appeared  was  determined  to  lie  between  137  and  168  m/s.  Unfortu¬ 
nately  the  rod  was  yawed  in  shots  103  and  104,  and  the  negative 
in  shot  105  turned  out  to  be  unaccountably  light  so  that  no 
pictures  are  available  that  show  complete  rod  profiles  for  normal 
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TABLE  8 


Shot  # 
0094 

0095 

0096 

0097 

0098 

0099 

00100 

0101 

0102 

0103 

0104 

0105 


TAYLOR  TEST  SHOT  MATRIX 


Velocity 

m/s 


Comments 


344 

293 

225 

169 

197 

197 

198 

217 

215 

138 

158 

133 


Good  impact,  radial  cracks, 
tip  diameter  =  0.710" 

Good  impact,  radial  cracks, 
tip  diameter  =  0.620" 

Good  impact,  radial  cracks, 
tip  diameter  =  0.687" 

Projectile  yawed,  no  radial 
cracks,  frame  rate  =  500,180  Hz 

Projectile  yawed,  one  radial 
crack,  frame  rate  =  500,100  Hz 

Good  impact,  no  radial  cracks, 
tip  diameter  =  0.450" 

Good  impact,  no  radial  cracks, 
tip  diameter  =  0.454" 
frame  rate  =  500,000  Hz 

Good  impact,  no  radial  cracks, 
tip  diameter  =  0.455",  frame 
rate  =  5000,000  Hz 

Good  impact,  one  radial  crack, 
tip  diameter  =  0.472" 

Good  impact,  no  cracks, 
tip  diameter  =  0.380", 
pictures  underexposed. 

Good  impact ,  no  cracks 
tip  diameter  =  0.405", 
projectile  yawed. 

Projectile  yawed. 


« 
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Figure  78. 


Post-impact  Photographs  of  Armco  Iron  Rods. 
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impacts.  At  this  stage,  the  experiments  were  discontinued 
because  no  more  Armco  iron  samples  were  available. 

The  images  from  a  total  of  four  shots:  0094,  0100, 
and  0101,  were  digitized  and  the  diameters  of  the  rods  as  a 
function  of  time  computed.  As  the  rod  was  crushed  during  the 
impact,  the  number  of  diameter  measurements  naturally  decreased. 
All  length  measurements  are  relative  to  the  back  end  of  the  rod. 
The  number  of  places  where  diameter  measurements  could  be  made 
also  diminishes  as  the  rod  advances  due  to  the  rising  debris 
cloud.  Figure  80  presents  the  data  for  these  measurements. 
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Figure  80.  Profiles  Determined  by  Digitizing  Successive 
Frames  in  Shot  94. 


SECTION  3 

CONSTITUTIVE  MODELING 


The  inatheinat ical  description  of  the  relationship  between 
stresses,  strains,  and  their  time  derivatives  is  referred  to  as 
the  constitutive  relationship  of  the  material.  Correct  and 
complete  constitutive  models  describe  material  behavior  over  the 
entire  range  of  states  that  may  be  experienced  in  impact  and 
explosive  loading.  However,  this  may  be  extremely  difficult  to 
achieve  (even  for  uniaxial  stress)  and  therefore  most 
constitutive  relationships  are  applicable  only  for  a  narrow  range 
of  stress  and  strain-rate  conditions.  The  most  elementary 
constitutive  models  are  appropriate  only  for  a  rate-independent 
material  under  uniaxial  monotonic  loading  conditions.  Perzyna 
(Reference  39),  Malvern  (Reference  40),  and  several  other 
investigators  proposed  various  visco-plastic  constitutive  models 
for  describing  the  material  response  at  high  strain-rate  loading 
conditions  [>1000  s  ^] .  In  those  models,  stress  can  be  computed 
for  a  given  strain-rate  and  strain,  but  load  history  dependency 
of  the  material  response  can  not  be  explicitly  handled. 

Several  constitutive  theories  based  on  state  variables  have 
been  developed  by  Bodner  and  Partom  (Reference  4),  Krieg  et  al., 
(Reference  41),  and  Chaboche  (Reference  42).  One  of  the 
advantages  of  state  variable  theories  is  that  changes  in 
constitutive  behavior  with  load  history  can  be  it.O(  eled  oy 
suitably  modifying  the  evolution  equation  of  the  corresponding 
state  variable.  Depending  on  the  complexities  of  the  loading, 
more  than  one  state  variable  may  be  introduced  into  tne  tnecry. 

For  the  present  work,  the  state  variable  based  const ivucive 
theory  developed  by  Bodner  and  Partom  was  selected  for  three 
reasons:  (1)  ability  to  predict  the  response  of  material  to  a 

.  road  range  of  load  histories;  (2)  aoaptability  to  finite 
difference  analyses  of  structural  components;  and  (3)  ease  with 
w'lich  the  parameters  in'  the  constitutive  models  can  be  determined 
from  high  strain-rate  test  data.  The  model  has  been  incorporated 


in  the  MAGNA  finite  element  code  and  the  STEALTH  finite 
difference  code. 


3.1  BODNER-PARTOM  MODEL 


For  assumed  small  strains,  the  total  strain  rate  is  taken 
to  be  decomposable  into  elastic  (reversible)  and  inelastic  (non- 
reversible)  components 


•e 

=  e  .  . 
ij 


+  eP 


which  are  both  non-2ero  for  all  loading/unloading  conditions. 
Inelastic  stresses  and  strains  corresponding  to  time  dependent 
reversible  deformations  with  energy  losses  (anelastic  behavior) 
are  not  treated  in  this  formulation  because  they  are  considered 
to  be  relatively  unimportant.  They  could  be  introduced  at  a 
later  stage  to  account  for  certain  transient  effects  associated 
with  creep  and  stress  relaxation. 

The  elastic  strain  rate,  j f  is  related  to  the  stress  rate 
by  the  time  deviation  of  Hooke's  Law  (based  on  the  assumption  of 
small  strains).  The  inelastic  strain  rate  e . .  can  be  expressed 

1 3 

in  the  general  form 

eP.  =  eP.(o.  .,  Z.  ,T)  (23) 

IJ  IJIJK 

where  Z  are  one  or  more  internal  (inelastic)  state  variables  and 
T  is  the  temperature.  In  particular,  is  taken  to  follow  the 

Prandtl-Reuss  flow  law  of  classical  plasticity,  (assuming 
incompressibility) : 

eP.  =  e?.  =  XS..  (24) 

11  11 

where  e*;  .  and  S..  are  the  deviatoric  plastic  strain  rate  and  the 
11  11 

stress.  Squaring  equation  (24)  gives 


1/2  eP.  e?.  =  dP  =  1/2\^S .  .  S. .  =  X^J 
11  11  2  '  11  11  2 
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(25) 


^2  =  ^^^2 


where  is  the  second  invariant  of  the  plastic  strain  rate  and 
is  the  second  invariant  of  the  stress  deviator. 

The  basic  assumption  of  the  formulation  is  that  all 
inelastic  deformations  are  governed  by  a  continuous  relation 
between  and  J^.  The  particular  form  that  was  established  to 
describe  this  relation  was  based  on  the  concepts  and  equations  of 
dislocation  dynamics  and  is  given  by: 

r^P  r^2  r  /  ^  (n  +  l){  /0/-X 

=  Do  e=<p(--35-)  —  (26) 


where  is  the  limiting  value  of  the  plastic  strain-rate  in 
shear.  n  is  a  material  constant  that  controls  strain  rate 
sensitivity  and  also  influences  the  overall  level  of  the  stress- 
strain  curves.  It  is  a  fundamental  constant  and  it  is  not 
dependent  on  the  loading  history.  Z  is  the  inelastic  state 
variable  which  represents  the  measure  of  the  overall  resistance 
to  plastic  flow  caused  by  microstructural  barriers  that  impede 
dislocation  motion.  Combining  equations  (24)  through  (26),  the 
plastic  strain-rate  takes  the  form 


e?  .  =  D  exp 
1  j  o  ^ 


(t?  V  (n-H)  ' 

~\^^2)  2n  J 


(27) 


The  general  form  for  the  evolution  equations,  i.e.,  history 
dependence,  of  the  inelastic  state  variable  Z  is, 

Z  =  F(J2,  Z,  T)  (28) 


For  conditions  of  uniaxial  stress  of  constant  sign,  the 
strain  hardened  material  is  assumed  to  be  represented  by  a  single 
state  variable  Z  which  depends  on  plastic  work.  This  corresponds 
to  isotropic  hardening  and  the  evolution  equations  employed  can 
be  shown  to  be  in  the  form  of  equation  (28).  The  corresponding 
form  is  then  based  on  the  concept  that  plastic  work,  controls 

the  hardening  process,  and  that  the  plastic  work  and  its  time 
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derivative  are  functions  of  a.,  and  e...  The  evolution  equation 

13  ij 

for  the  state  variable  is  assumed  to  have  the  form: 

Z  =  m(Z, -Z)  W  (29) 

1  P 

where,  W  =0..  . .  The  above  equation  can  be  integrated  with  m 

P  ^  J  ^1 

and  Z,  as  material  constants.  The  relation  between  Z  and  W 

P 

becomes , 

Z  =  Z^-(Z^-Z^)  e""*^p  (30) 

Z^  is  the  initial  value  of  Z  corresponding  to  the  reference  state 

from  which  W  is  measured.  It  is  noted  that  Z  could  take  any 
p  o 

value  between  0  and  The  initial  yield  stress  of  the  uniaxial 

stress-strain  relation  is  directly  influenced  by  Z^.  Z^ 
corresponds  to  the  saturation  state  of  Z  at  which  the  material 
reaches  a  fully-work  hardened  condition.  m  is  a  material 
constant  that  controls  the  rate  of  work  hardening. 

For  materials,  such  as  OFHC  copper,  which  exhibit  strong 
strain  hardening,  m  was  made  a  function  of  by  Bodner  and 
Partom  (Reference  4).  The  expression  for  m  was: 


m  =  m  +  m,  e 
o  1 


which  adds  two  more  constants.  When  m  is  defined  by  equation 

(31),  the  relationship  (30)  is  not  valid.  The  expression  for  Z 

when  integrated  after  replacing  m  by  the  definition  (31),  yields 

the  following  relationship  between  Z  and  W  : 

P 

„  -(m  +  m^  -  m) 

-m  W  _o _ 1 _ 

Z  =  Zi-(Zi-  Zq)  e  °  P*e  a  (32) 


For  evaluating  the  BP-model  constants  from  the  SHB  and 
plate  impact  test  results,  we  need  to  develop  1)  a  computer 
program  to  solve  the  BP-model  equations  under  uniaxial  stress 
condition  (SHB  test),  and  2)  special  purpose  subroutines  for 
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STEALTH  to  describe  the  BP-inodel  under  one  dimensional  strain 
condition  (plate-impact  test). 

A  program  called,  BPSOLVE  was  developed  to  calculate  the 
uniaxial  stress  for  a  given  strain-rate  history.  This  program 
was  effectively  used  in  the  BP-model  constants  evaluation.  The 
following  section  describes  the  efforts  that  were  taken  to 
program  the  BP-model  in  STEALTH. 


3.2  BODNER-PARTOM  MODEL  IN  STEALTH 

Currently,  the  standard  STEALTH  code  is  not  capable  of 
describing  a  visco-plastic  material.  Since  the  state  variable 
based  Bodner-Partom  constitutive  equations  were  selected  for  high 
strain-rate  material  characterization,  it  was  essential  to 
incorporate  the  BP-model  algorithm  into  the  STEALTH  code  for  the 
plate  impact  test  simulations  using  the  BP-model  parameters.  To 
incorporate  the  BP-model  into  STEALTH,  various  subroutines  were 
developed . 

A  subincremental  time  stepping  scheme  was  used  to  calculate 
the  deviatoric  stresses,  S^^  for  each  finite  difference  time 
step.  The  global  time  step.  At  is  further  divided  into  small 
steps  in  a  special  purpose  subroutine  to  describe  the  BP-model. 

At  the  end  of  each  global  time  step,  STEALTH  provides  the 

•  •  p 

volumetric  strain  rate,  e,  and  the  pressure,  p;  S. o. e . .  are 

t'  '  '  2.  j  ij  ij 

unknown . 


In  the  current  numerical  scheme,  the  plastic  strain-rates 
were  estimated  for  each  subincremental  time  step  using  equation 
(25)  with  the  value  for  J„  calculated  from  the  known  S. .  values 

2  ij 

•  © 

at  the  beginning  of  the  step.  The. elastic  strain-rates,  e. ., 

*  D  ^ 

were  estimated  from  the  estimated  £5.  and  the  known  total  strain 
.  .13 

rate,  e. ..  The  new  estimates  for  S..  can  be  made  from  the 
13  13 

following  relationship: 


where 


S.  . 

13 


•  • 

a .  .  +  p 

13 


a .  . 

13 


•  e 

e  .  . 


13 


(33) 

(34) 
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and 


E. .  are  the  elements  of  the  elastic  modulus 

matrix. 


Using  the  new  estimates  for  j  f  the  improved  values  for 
the  plastic  strain-rates  at  the  end  of  the  subincremental  time 
step  can  be  computed.  This  procedure  is  continued  until  the 
values  of  computed  for  the  two  successive  iterations  converge 

to  the  same  value  within  the  tolerant  limits.  A  series  of 
flowcharts,  given  in  Appendix  A  describe  the  steps  involved  in 
the  numerical  scheme. 


Rajendran  and  Grove  (Reference  43)  discuss  in  detail  the 
various  newly  developed  subroutines  for  STEALTH  and  also  the 
methods  that  were  used  for  the  validation  of  the  special  purpose 
routines . 


4 


9, 


3.3  EVALUATION  OE  BODNER-PARTOM  MODEL  PARAMETERS 

The  Bodner-Partom  constitutive  model  contains  seven 

empirically  derived  constants;  D^,  n,  m^,  m^,  a,  and  Z^;  these 

constants  were  defined  in  Section  3.1.  D  is  a  reference  maximum 

o 

strain-rate  level  which  can  be  conveniently  assigned  a  large 

value.  In  this  investigation,  D  was  arbitrarily  chosen  to  be 
8—1  ^ 

10  s  .  The  constants  n,  Z^,  and  m(=m^)  were  determined  by 

first  estimating  the  values  based  on  the  Hopkinson  bar  data,  then 
modifying  the  values  so  that  the  stress-strain  predictions  from 
the  model  match  both  the  Hopkinson  bar  and  the  plate  impact  test 
results  as  closely  as  possible.  Eor  most  materials,  m  was 
assumed  to  be  a  constant  and  represented  by  strong 

strain-hardening  materials  like  OFHC  copper,  m^^,  and  a  had  to  be 
added  to  the  other  material  constants. 

•  n 

The  plastic  strain  rate,  under  uniaxial  stress  condition 
is  given  by 
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The  plastic  strain  rate  (approximately  equal  to  the  total 

•  0  •  n 

strain  rate,  since  e  <<  e  )  and  the  maximum  stress  (at  which  the 
stress-strain  curve  is  approximately  flat)  can  be  obtained  from 
each  Hopkinson  bar  tensile  test.  At ,  a  =  a  ,  the  state 
variable  Z  saturates  and  reaches  its  maximum  value  (Z  =  Z^). 

Then,  equation  (35)  can  be  rewritten,  after  taking  natural 
logarithms  twice,  as; 


where 


and 


Y  =  M  X  +  C 

(36) 

/3  eP 

Y  =  ln[-ln(---^3-)] 
o 

(37) 

X  =  ln(  a  ) 

max 

(38) 

c 

CN 

1 

II 

(39) 

C  =  ln((  2^1  q  1 

(40) 

Since  various  values  of  X  and  Y  are  known  from  different 
strain-rate  Hopkinson  bar  tensile  tests,  a  straight  line  can  be 
fitted  to  the  data  and  the  corresponding  constants  M  and  C  can  be 
obtained  from  the  slope  and  intercept  of  the  straight  line. 
Knowing  M  and  C,  the  Bodner-Partom  model  constants  n  and  Z^  can 
be  calculated  from  equations  (37)  and  (38). 

The  initial  estimates  of  the  constants  Z  and  m  were 

o 

determined  from  the  plot  of  plastic  work  (W  )  vs.  Z,  where  W  and 

P  P 

Z  are  defined  as; 

eP 

W  =  /  ad  eP  (41) 

P 

o 


and 


Z  =  [-ln( 


/3 


eP 


2  D 


)  ( 


n 


o 


n  +  1 


2ni l/2n 
a  ] 


(42) 
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To  allow  for  direct  integration  of  equation  (39),  the 
Hopkinson  bar  stress-strain  data  were  described  in  terms  of  a 
polynomial  curve  fit. 

0  =  f(eP)  (43) 

A  linear  fit  was  used  to  describe  the  plot  of  vs.  Z. 

W  =  A  Z  +  B  (44) 

P 

The  initial  estimates  for  m  and  Z^  are  determined  from  the  slope 
(A)  and  intercept  (B)  of  the  curve  fit. 

m  =  -A  (45) 


Since  regression  of  all  the  Hopkinson  bar  test  data  is  a 
time  consuming  and  a  tedious  process,  a  series  of  interactive 
automated  computer  programs  was  developed  to  evaluate  the 
material  constants.  A  flow  chart  for  the  interactive  program  is 
shown  in  Figure  81. 

Since  n  controls  the  flow  stress  level  due  to  strain  rate 

effects  and  Z  controls  the  initiation  of  plastic  flow,  these  two 
o 

parameters  greatly  influence  the  HEL  in  the  plate  impact  test 
simulations;  whereas  all  the  four  parameters  (m,  2^,  n  and  Z^) 
influence  the  stress  levels  in  the  SHB  test  simulations  at 
strains  >0.02.  Due  to  the  noise  in  the  SHB  test  data  at  small 
strains,  the  SHB  flow  stress  values  at  such  strains  are  not 
dependable.  Therefore,  it  is  essential  to  use  the  plate  impact 
test  data  (HEL)  in  the  evaluation  of  BP-model  parameters  at  small 
strains . 

Unfortunately,  in  the  interactive  program  it  is  not 
possible  to  include  the  flow  stress,  Y^,  obtained  from  the  HEL  of 
the  plate  impact  test  directly  because  the  strain-rate  associated 
with  the  HEL  is  not  a  measurable  quantity.  To  overcome  this 
difficulty  and  also  to  make  use  of  the  experimentally  obtained 


148 


81.  Flow  Chart  for  the  Interactive  Computer 
Program  Describing  Bodner-Partom  Model 
Constants  Evaluation. 


:e 


Figure  81  (cont'd).  Flow  Chart  for  the  Interactive  Computer 

Program  Describing  Bodner-Partom  Model 
Constants  Evaluation. 


Y^,  an  indirect  data  reduction  scheme  was  developed  and 
successfully  used  for  all  five  metals.  The  steps  involved  in 
that  scheme  were: 

1.  Estimation  of  n,  m  ,  Z  ,  Z,,  from  the 

o  o  1 

interactive  program  using  only  the  SHE 
tension  test  results. 

2.  Simulate  the  plate  impact  test  through  STEALTH 
using  the  newly  developed  BP-model  algorithm  to 
check  whether  the  HEL  based  on  the  estimated 
constants  predicted  the  test  data  correctly. 

If  the  predicted  HEL  is  above  the  test  result, 
the  value  for  n  must  be  adjusted  for  less  rate 
dependency  and  vice  versa . 

3.  Reduce  or  increase  the  value  of  n  according  to 
Step  2.  Simulate  the  SHE  test  results  using 
EPSOLVE  with  the  new  estimate  of  n.  Adjust 
other  parameters  interactively  until  SHE  test 
data  are  matched. 

4.  Perform  Steps  2  and  3  until  both  SHE  and  plate 
impact  test  data  agreed  with  the  model 

cal ibrat ions . 

Using  this  data  evaluation  scheme,  model  constants  for  the  five 
materials  were  successfully  obtained.  The  values  are  listed  in 
Table  9.  Following  sections  describe  the  model  predictions. 

3.3.1  Model  Evaluations 

This  section  shows  the  results  of  the  E-P  simulations 

of  tensile  SHE  and  plate  impact  tests.  We  also  demonstrate  the 

ability  of  the  model  to  predict  the  results  of  compression  SHE 

tests,  in  which  the  strain  rate  varies  continually.  For  most 

materials,  it  was  only  necessary  to  evaluate  the  constants  n,  m^ , 

Z  ,  and  Zt  of  the  Bodner-Partom  constitutive  model.  For  OFHC 
o  1 

copper  the  additional  constants  m^  and  a  were  needed  to  predict 

the  SHE  test  data  accurately.  The  constant  D  was  assumed  to  be 

8—1  ^ 

10  s  for  all  the  five  materials.  The  continuously  varying 
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TABLE  9 


B-P  MODEL  AND  ELASTIC  CONSTANTS 


1020 

OFHC 

C1008 

HYIOO 

7039-T64 

K  (Kbar) 

1540 

1418 

1540 

1540 

768 

G  (Kbar) 

794 

450 

794 

794 

276 

E  (Kbar) 

2033 

1220 

2033 

2033 

739 

V 

0.28 

0.36 

0.28 

0.28 

0.34 

(Kbar) 

4.6 

1.6 

6.9 

9.5 

3.7 

Do  (s'L 

10« 

10® 

10® 

10® 

10® 

n 

4.0 

0.4 

0.4 

1.2 

4.0 

in  (Kbar  ^) 
o 

3.0 

1.1 

1.5 

1.0 

2.8 

(Kbar  ^) 

0 

15.0 

0 

0 

0 

Z  (Kbar) 
o 

6.4 

8.0 

55. 

24. 

5.6 

(Kbar) 

9.3 

65.5 

70. 

35.5 

7.6 

a  (Kbar  ^) 

0 

150 

0 

0 

0 

152 


strain-rate  data  obtained  from  the  tests  were  used  as  program 
input  in  the  compression  test  simulations. 

a .  1020  steel 

The  Bodner-Partom  model  constants  for  this 

material  were  obtained  based  on  the  selected  SHB  test  results. 

It  was  observed  that  that  Y  (=4.6  Kbar)  obtained  from  the  plate 

o 

impact  tests  was  no  higher  than  the  yield  values  obtained  in  the 
SHB  tests.  This,  along  with  observations  summarized  by  Bless 
(Reference  19),  supported  the  assumption  that  1020  steel  is  a 
fairly  rate-independent  material.  Matuska  (Reference  44)  carried 
out  finite  difference  calculations  to  reproduce  hemispherical 
impact  data  of  Bless  and  Barber  (Reference  45)  for  this  material. 
He  had  to  employ  a  work  hardening  model  (with  =  4.6  Kbar  at 
initial  yield,  increasing  linearly  to  6.3  Kbar)  designed  to  fit 
the  observed  velocity  decay.  In  Figure  82  three  different  SHB 
tests  and  the  model  used  by  Matuska  are  shown.  It  clearly 
supports  the  conclusion  that  this  steel  is  a  fairly  rate 
independent  material  in  the  high  strain-rate  regime. 

The  constants  of  the  Bodner-Partom 
constitutive  model  for  the  true  stress-true  strain  behavior  of 
1020  steel  are  summarized  in  Table  9.  Figure  83  contains  the 
actual  test  results  and  the  Bodner-Partom  predictions  for  each 
strain  rate  level  tested. 

The  constants  of  the  Bodner-Partom 
constitutive  model,  derived  from  the  Hopkinson  bar  tensile 
experiment  and  the  plate  impact  tests,  were  used  to  predict  the 
dynamic  compressive  true  stress-true  strain  behavior.  The 
predictions  and  the  original  Hopkinson  bar  data  are  shown  in 
Figure  84.  The  slight  disagreement  shows  that  there  is  some 
anisotropy  present  in  the  materials.  However,  the  BP-model 
precision  is  adequate  for  engineering  calculations. 

Unfortunately,  none  of  the  VISAR  plots  for 
this  material  turned  out  to  be  useful.  Among  the  various  shots, 
shot  #529  happened  to  be  slightly  better  than  the  rest. 

Simulation  of  this  test  was  attempted  through  STEALTH  with  the 


Figure  82. 


Three  Different  SHE  Tests  and  the  Model 
Used  by  Matuska  (Reference  44)  . 
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Bodner-Partoin  Prediction 

A/i 

Hopkinson  Bar  Results 


TRUE  STRAIN 
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Hopkinson  Bar  Results 
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d.  Combined  Plot 


Figure  83.  Bodner-Partom  Predictions  of  the  Dynamic  Tensile 
Behavior  of  1020  Steel. 


Bodner-PartOTTi  model  constants  in  Table  10  for  1020  steel.  The 
VISAR  trace  of  the  free  surface  velocity  history  was  compared  with 
the  simulation  in  Figure  85.  Since  the  shape  of  the  rising 
portion  of  the  VISAR  beyond  the  HEL  was  so  poor,  it  is  difficult 
to  judge  the  accuracy  of  model  predictions. 

b.  QFHC  copper 


The  constants  of  the  Bodner-Partom 
constitutive  model  for  the  true  stress-true  strain  behavior  of 
OFHC  copper  are  summarized  in  Table  9.  Figure  86  contains  the 
actual  SHB  test  results  and  the  Bodner-Partom  simulation  for 
several  strain  rate  levels.  Since  this  material  exhibited 
substantial  strain  hardening,  modeling  with  only  four  constants, 
n,  m^,  Z^,  and  was  not  successful.  As  Bodner  and  Partom 
(Reference  4)  reported  earlier,  additional  constants  m^ ,  and  a 
were  reguired  in  the  modeling.  The  plate  impact  test  results  were 
not  used  in  the  model  calibration  because  of  the  uncertainty  in 
the  Since,  the  HEL  value  was  so  low,  it  was  not  well 

resolved.  However,  a  value  of  1.6  Kbar  reported  earlier 
(Reference  19)  seems  to  support  the  model  predictions. 


The  constants  of  the  Bodner-Partom 
constitutive  model,  derived  from  the  Hopkinson  bar  tensile 
experiments,  were  used  to  predict  the  dynamic  compressive  true 
stress-true  strain  behavior.  The  predictions  and  the  original 
Hopkinson  bar  data  are  shown  in  Figure  87.  The  disagreement 
indicates  that  the  OFHC  copper  exhibits  anisotropic  behavior. 


Excellent  VISAR  records  for  OFHC  copper  plate 
impact  experiments  are  given  in  Bless  (Reference  19)  and  Bless  and 
Paisley  (Reference  46).  For  model  prediction,  shot  #560  (4mm 
flyer  and  6mm  target  at  209  m/s)  was  simulated  through  STEALTH 
using  the  BP-model  subroutines  to  describe  the  high  strain-rate 
material  behavior  and  the  constants  are  given  in  Table  9.  The 
model  prediction  of  the  free  surface  velocity  of  the  target  is 
compared  directly  with  the  VISAR  plot  in  Figure  88.  Since  the  HEL 
of  copper  is  very  low,  it  is  difficult  to  compare  the  HEL;  however, 
the  plots  beyond  the  HEL  which  includes  the  peak  velocity  and  the 
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Figure  85.  VISAR  Trace  and  Bodner-Partoin  Prediction 
for  1020  Steel. 
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Figure  86.  Bodner-Partom  Prediction  of  the  Dynamic  Tensile 
Behavior  of  OFHC  Copper. 


Figure  86  (cont'd).  Bodner-Partom  Predictions  of  the  Dynamic 

Tensile  Behavior  of  OFHC  Copper. 


VELOCITY  (KM/SEC) 


shape  of  the  release  portion  compared  extremely  well.  A  typical 
STEALTH  plot  of  the  stress-strain  history  in  the  spall  plane  is 
shown  in  Figure  89.  The  complex  dynamic  stress  history  at  the 
spall  plane  can  be  calculated  with  great  detail  by  STEALTH  due  to 
the  sophisticated  constitutive  model  (BP-model)  description. 

The  similar  SWAP  code  was  also  used  for  an 
analysis  of  OFHC  copper.  The  objective  was  to  infer  the  uniaxial 
stress-strain  curve  that  is  consistent  with  the  release  data  from 
shot  560.  This  code  can  include  strain  hardening,  but  not  strain 
rate  dependency.  The  approach  was  to  incrementally  vary  the 
hardening  model  while  continually  comparing  predicted  and  actual 
release  wave  profiles.  It  was  found  that  the  release  wave  profile 
in  these  relatively  thick  targets  was  very  sensitive  to  the 
assumed  hardening  model. 

Hardening  is  introduced  into  SWAP  by  making 
the  yield  strength  a  function  of  the  plastic  strain.  The 
equation  used  was 

Y  =  +  ae^  (4  7) 

o 

where  a  and  b  are  empirical  constants. 

The  effect  of  changing  the  hardening  behavior 
is  shown  in  Figures  90  and  91.  Both  figures  show  the  surface 
velocity  profiles  for  shot  560.  In  Figure  90,  the  material  is 
assumed  to  be  elast ic-perf ect ly  plastic  with  a  yield  strength  of 
1.0  Kbar  (a  =  0).  In  Figure  91,  the  material  is  assumed  to  yield 
at  the  same  initial  stress  with  a  =  3.4  Kbar  and  b  =  0.28. 

The  figures  include  a  comparison  with  the  observed  velocity 
profile.  The  agreement  is  excellent,  until  spall  occurs.  (Spall 
was  suppressed  in  the  calculation. )  Similar  good  results  were 
obtained  for  other  shots  with  this  model. 

Figure  92  presents  the  uniaxial  stress-strain 
curve  implied  by  these  parameters.  This  curve  represents  the 
average  stress  strain  relationship  pertaining  to  a  tensile  test 
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Figure  89.  The  Effective  Stress  vs.  Effective  Strain 
at  the  Spall  Plane  for  OFHC  Copper  from 
STEALTH  with  Bodner-Partom  Routines. 


Figure  90.  Comparison  of  SWAP  with  VISAR  Data  for  OFHC 
Copper,  No  Hardening. 
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Figure  91.  Comparison  of  SWAP  with  VISAR  Data  for  OFHC 
Copper,  Using  Hardening  Shown  in  Figure  92. 
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Figure  92.  Uniaxial  Stress-Strain  Curve  for  OFHC  Copper 
Derived  from  SWAP  Simulation. 
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at  ultra  high  strain  rates.  It  lies  considerably  above  even  the 
highest  strain  rate  SHB  test  data  (see  Figure  86e). 


c.  C1008  steel 

The  constants  of  the  Bodner-Partom 
constitutive  model  for  the  true  stress-true  strain  behavior  of 
C1008  steel  are  summarized  in  Table  9.  Figure  93  contains  the 
actual  test  and  the  Bodner-Partom  model  results  for  each  strain 
rate  level  tested. 

The  constants  of  the  Bodner-Partom 
constitutive  model,  derived  from  the  Hopkinson  bar  tensile 
experiments,  were  used  to  predict  the  dynamic  compressive  true 
stress-true  strain  behavior.  The  predictions  and  the  original 
Hopkinson  bar  data  are  shown  in  Figure  94.  The  slight  disa¬ 
greement  shows  that  there  is  some  anisotropy  present  in  the 
materials.  However,  the  BP-model  precision  is  adequate  for 
engineering  calculations. 

The  plate  impact  test  of  shot  #678  was 
simulated  through  STEALTH.  A  3mm,  1020  steel  flyer  and  6mm, 

C1008  steel  target  were  modeled.  The  impact  velocity  was  178 
m/s.  The  simulated  free  surface  velocity  history  is  compared 
with  the  VISAR  plot  in  Figure  95.  The  observed  HEL  is  slightly 
higher  than  the  calculated  value.  This  is  consistent  with  the 
precursor  decay  discussed  in  Section  2.3.3. 

The  STEALTH  simulation  reproduces  the  ramp 
on  loading  and  the  ramp  behind  the  plastic  wave  arrival.  This 
latter  feature  is  apparently  caused  by  rate  dependency.  However, 
the  peak  velocity  is  not  compared  well.  To  check  this  point, 
particle  velocity  ('-'p)  vs.  stress  [a)  diagram  was  constructed 
using  the  shock  relationship: 


a 


%EL 


^HEL 


(48) 


where 


^HEL' 


particle  velocity 


and  u  are  the  stress,  density  and 

Po 

at  the  elastic  limit,  respectively. 


the 

The 


Figure  93.  Bodner-Partoin  Predictions  and  the  Original 
Hopkinson  Bar  Results  for  Tensile  SHB  Tests 
of  C1008  Steel. 
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Figure  93  (cont'd) .  Bodner-Partom  Predictions  and  the  Original 

Hopkinson  Bar  Results  for  Tensile  SHB  Tests 
of  C1008  Steel. 


c.  High  Strain  Rate  Range 


Figure  94  (cont'd).  Bodner-Partom  Predictions  of  the  Dynamic 

Compressive  Behavior  of  C1008  Steel. 
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Figure  95.  Comparison  of  Velocity  History  Between 
STEALTH  with  Bodner-Partom  Routines  and 
VISAR  Data  (C1008  Steel). 
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diagram  constructed  for  shot  #678  is  given  in  Figure  96.  For  the 
impact  velocity  of  178  m/s,  the  calculated  peak  velocity  of  the 
target's  free  surface  was  0.1 63mm/Msec.  It  can  be  seen  from 
Figure  95  that  the  velocity  in  the  simulation  ramps  to  the  same 
peak  velocity.  The  peak  velocity  of  the  VISAR  plot  was  around 
0.174mm/ysec.  The  difference  between  the  theoretically  obtained 
peak  velocity  and  the  VISAR  can  be  attributed  to  uncertainty  in 
the  VISAR  data  reduction.  Due  to  the  discrepancy  in  the  peak 
velocities,  the  release  characteristics  of  the  plots  in  Figure  95 
did  not  match  well;  however,  the  slopes  and  the  trends  are  very 
similar.  It  appears  from  these  plots,  that  the  VISAR  data 
corrected  to  the  theoretically  predicted  peak  velocity,  may 
compare  well  with  the  simulation. 

To  validate  the  BP-model  constants  further, 
shot  #  685  (u^  =  134  m/s)  was  simulated  and  compared  with  the 
VISAR  plot  in  Figure  97.  The  results  were  similar.  The  HEL  was 
slightly  underprecicted  due  to  precursor  decay.  The  ramped  shock 
wave  was  reproduced.  Additional  experiments  should  be  conducted 
with  thin  material  to  verify  the  shape  of  the  ramped  wave.  This 
is  a  usual  practice.  Tests  are  needed  to  verify  its  reality  and 
determine  the  impact  conditions  that  give  rise  to  its 
occurrence. 

d.  HYIOO  steel 

The  constants  of  the  Bodner-Partom 
constitutive  model  for  the  true  stress-true  strain  behavior  of 
HYIOO  steel  are  summarized  in  Table  9.  Figure  98  contains  the 
original  Hopkinson  bar  results  and  the  Bodner-Partom  predictions 
for  each  strain  rate  level  tested. 

The  constants  of  the  Bodner-Partom 
constitutive  model,  derived  from  the  Hopkinson  bar  tensile 
experiments,  were  used  to  predict  the  dynamic  compressive  true 
stress-true  strain  behavior.  The  predictions  and  the  original 
Hopkinson  bar  data  are  shown  in  Figure  99.  Rate  dependency  in  the 
model  is  slight,  and  the  agreement  with  predicted  and  measured 
flow  stress  is  within  10  percent. 


STRESS  (Kbars) 


Figure  96.  Stress  vs.  Particle  Velocity  Diagram  for 
C1008  Target  and  1020  Flyer  Plate. 
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Figure  97 


Comparison  of  Velocity  History  Between  STEALTH 
with  Bodner-Partom  Routines  and  VISAR  Data  for 
C1008  Steel  (Impact  Velocity  =  134  m/s) . 
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Figure  98. 


Bodner-Partom  Predictions  and  the  Original 
Hopkinson  Bar  Results  for  Tensile  SHB  Tests 
of  HYIOO  Steel. 
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The  plate  impact  test  of  shot  #680  was 
simulated  with  STEALTH  with  the  same  BP-model  constants  that  were 
used  in  the  SHB-test  simulations.  The  impact  velocity  was  355 
m/s.  The  simulated  free  surface  velocity  plot  compared  reasonably 
well  with  the  VISAR  plot  as  shown  in  Figure  100.  The  HEL  and  the 
release  characteristics  of  these  two  plots  match  extremely  well. 
The  peak  velocity  in  the  simulation  is  within  10%  of  the  observed 
value.  Since  theoretical  estimation  of  the  peak  particle  velocity 
using  the  equation  (48)  coincided  with  the  velocity  simulation, 
one  believes  that  the  difference  may  be  caused  by  drift  in  the 
VISAR  signal  (e.g.,  slow  change  in  fringe  contrast). 

e .  7039-T64  aluminum 

The  Bodner-Partom  constants  for  7039-T64 
aluminum  are  summarized  in  Table  9.  Figure  101  contains  the 
original  Hopkinson  bar  results  and  the  Bodner-Partom  simulations 
for  each  strain  rate  level  tested.  Since  the  material  is 
relatively  rate  independent,  n  takes  on  a  large  value.  Based  on 
the  tensile  results,  the  Bodner-Partom  model  was  used  to  predict 
results  of  compressive  true  stress-true  strain  measurements.  The 
predictions  and  the  original  Hopkinson  bar  data  are  shown  in 
Figure  102.  The  predictions  are  generally  good  and  indicate 
little  anisotropy  in  the  material. 

Shot  #682  was  simulated  using  the  BP-model 
constants  given  in  Table  9  for  7039-T64.  The  impact  velocity  was 
375  m/s.  The  simulation  is  compared  with  the  VISAR  plot  in  Figure 
103.  The  HEL,  the  rising  portion  beyond  HEL,  the  peak  velocity 
and  also  the  release  characteristics  of  the  simulation  compared 
extremely  well  with  the  VISAR.  The  accuracy  with  which  both  the 
SHB  and  plate  impact  test  were  simulated  strengthens  confidence 
in  the  values  for  the  B-P  constants  and  indicates  the 
dependability  of  the  modeling  techniques. 

3.4  SPALL  FAILURE  PARAMETERS 

Spall  damage  is  also  a  function  of  the  amplitude  of  the 
tensile  stress  and  the  stress  duration.  Often  experiments  are 
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Figure  100.  Comparison  of  Velocity  History  Between 
STEALTH  with  Bodner-Partora  Routines  and 
VISAR  Data  for  HYlOO  Steel. 
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Figure  101.  Bodner-Partom  Model  Predictions  and  the  Hopkinson 

Bar  Results  for  Tensile  SHB  Tests  of  7039-T64  Aluminum. 
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c.  High  Strain  Rate  Range 


Figure  102. 


Bodner-Partom  Predictions  of  the  Dynamic 
Compressive  Behavior  of  7039-T64  Aluminum. 
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Figure  103.  Comparison  of  Velocity  History  Between  t 

STEALTH  with  Bodner-Partom  Routines  and 
VISAR  Data  for  7039-T64  Aluminum. 
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performed  with  various  thickness  flyer  plates  so  that  the 
amplitude/duration  criterion  for  spall  may  be  explicitly 
developed. 

Thus,  there  are  several  different  ways  to  specify  spall 
criteria: 

>  1.  The  spall  stress,  a  ,  given  by  equation  (13).  This 

s 

stress  probably  represents  the  physical  stress  in  the 
r  target  associated  with  initial  void  formation.  It  is 

normally  determined  from  experiments  in  which  the 
spall  is  complete. 

2.  The  peak  tensile  stress  in  a  shot  that  causes 
incipient  spall  failure.  This  is  a  function  of  the 
tensile  pulse  duration.  It  is  plausible  that  for 
thick  enough  flyer  plates,  this  stress  would  be  equal 
to  ,  although  that  has  not  been  experimentally 
verified.  For  steels,  there  is  relatively  little 
change  in  this  stress  for  flyer  plates  3mm  thick  or 
greater . 

3.  The  spall  parameter  that  predicts  spall  damage  in 
elementary  numerical  models  which  treat  spall 
fracture  as  instantaneous.  When  the  spall  threshold 
is  greatly  exceeded,  the  numerical  models  become  more 
accvirate,  since  spall  is  very  rapid.  However,  for 
near  threshold  conditions,  the  model  parameters  do 
not  correspond  to  physical  stresses,  since  the  actual 
mechanics  of  spall  damage  (i.e.,  gradual  void  growth 
and  coalescence)  is  not  described  by  the  model.  The 
spall  models  used  in  SWAP,  STEALTH,  and  most 

*  practical  finite  difference  codes  are  in  this 

category. 

Table  1  lists  several  spall  criteria.  The  first  four 
columns  give  parameters  that  correctly  describe  spall  when  used 
with  simplified  material  models.  Those  values  are  recommended  for 
numerical  simulations.  The  last  column  is  the  estimate  of  the 
most  tensile  stress  that  the  material  can  actually  sustain.  These 
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are  the  only  values  that  represent  material  properties.  They 
should  be  used  in  developing  theories  for  void  nucleation,  growth, 
and  coalescence. 

Present  efforts  were  not  focussed  on  modeling  the  spall 
failure  through  application  of  the  Bodner-Partom  constitutive 
model.  Instead,  for  engineering  applications,  two  spall  criteria  * 

were  selected  to  simulate  a  plate  impact  test  in  which  the  target 
plate  was  partially  spalled.  A  simple  material  model  based  on  , 

elast ic-perfectly  plastic  material  behavior  was  used  with  the 
experimentally  obtained  for  describing  the  flyer  and  target 
plates.  The  first  criterion  was  time  independent  and  based  on  a 
critical  spall  stress.  When  the  tensile  stress  reaches  a  critical 
level  ,  the  spall  fracture  occurs  instantaneously.  The  second 
criterion  was  time  dependent  and  based  on  a  critical  level  of 
damage.  The  Tuler  and  Butcher  (Reference  47)  formulation  was 
employed.  The  criterion  was  given  by; 

=  /  (0  -  0^)^  dt  (49) 

2,  ^xc  spall  (50) 

where  K^,  0^  and  X  are  fracture  parameters.  According  to  this 

model,  there  is  a  minimum  normal  stress,  0^,  at  which  spall  can 

occur.  At  0=  0^,  the  tensile  stress  must  be  applied  for  an 

infinitely  long  time  for  a  spall  plane  to  form.  When  0>  spall 

can  occur  even  for  a  short  pulse  duration,  depending  on  the  impact 

velocity  and  the  thickness  of  the  flyer  plate.  For  large  values 

of  ( 0  -  0^)  spall  takes  place  very  rapidly  and  the  time  dependence 

is  less  manifest.  Under  such  circumstances  it  often  turns  out 

that  spall  can  be  described  by  a  single  spall  stress  0^,  thus,  the  * 

first  criterion  may  be  sufficient,  according  to  Zukas  (Reference  1), 

Bertolf  et  al. ,  (Reference  48),  and  Bless  and  Paisley  (Reference  46).  * 

The  Tuler-Butcher  model  was  programmed  into  both  STEALTH  and  SWAP 

as  a  user  supplied  failure  criterion. 
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a 


1020  steel 


Bless  (Reference  19)  has  presented  a 
discussion  on  the  spall  fracture  of  this  steel.  In  his  work, 
spall  parameters  were  evaluated  through  SWAP  (Reference  49)  and 
HULL-code  simulations.  The  present  work  used  the  value  for 
(=4.6  Kbar)  as  reported  by  Bless  (Reference  19).  The  shot  #529 
was  simulated  through  STEALTH.  The  critical  spall  stress  and  the 
Tuler-Butcher  criteria  were  used  in  the  evaluations  of  the  spall 
parameters.  The  simulations  and  the  VISAR  data  are  compared  in 
Figure  104.  Figure  105  shows  the  rapid  fall  of  stress  preceding 
fracture  on  the  spall  plane.  The  value  of  2  was  assumed  for  the  A 
as  reported  by  Bless.  The  threshold  stress  in  the  Tuler-Butcher 
model  was  assumed  to  be  the  o  .  The  K  was  calculated  so  that 

*  ri  Ej  Li  A 

the  arrival  time  of  the  spall  signal  in  the  simulation  matches  the 
experiment.  K  =  14  Kbar  -  ys  was  found  to  be  the  best  value. 

A  more  complete  evaluation  of  the  Tuler- 

Butcher  parameters  was  conducted  with  the  SWAP  code.  Many 

calculations  were  performed  for  flyer  plate  thicknesses  between 

0.25  and  2mm.  Like  STEALTH,  SWAP  does  not  provide  for  growth  of 

damage.  Therefore,  it  is  inconsistent  to  use  this  code  to 

calculate  complete  spall,  since  the  material  properties  on  the 

spall  plane  are  affected  by  cracks  that  open  before  complete 

separation  takes  place.  Instead,  the  code  was  used  to  calculate 

an  incipient  spall  criterion.  Fracture  was  suppressed,  and  values 

of  K  were  computed  for  impacts  that  lie  on  the  no  spall/spall 

threshold.  The  best  agreement  between  K-contours  and  the  observed 

damace  threshold  was  for  A  =  2,  o  =  -10  Kbar,  and  K  ..  =  125 
2  '  o  crit 

Kbar  PS.  The  curve  for  these  values  is  drawn  in  Figure  56.  It 
was  calculated  with  d  =  2mm,  D  =  5mm;  d  =  1mm,  D  =  5mm;  and  d  = 
0.5mm,  D  =  2mm;  thereby  encompassing  a  slight  dependence  on  target 
thickness.  For  example,  at  d  =  1mm,  the  predicted  threshold  for 
D  =  2.5mm  is  about  20  m/s  lower  than  the  curve  as  drawn. 

The  Tuler-Butcher  parameters  used  in  Figure  56 
are  relatively  unique.  Smaller  values  of  result  in  curves  that 
rise  too  steeply  for  thin  flyer  plates,  while  larger  values  are 
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Figure  104.  Spall  Simulation  for  1020  Steel  Target  and 

Flyer  from  STEALTH  using  Tuler-Butcher  Spall 
Model  is  Compared  with  VISAR  Data. 
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Figure  105.  Stress  History  at  the  Spall  Plane  from 
STEALTH  Simulation  of  Shot  No.  529 
(1020  Steel) . 


inconsistent  with  shot  553.  Curves  with  X  =  3  also  rise  too 
steeply  as  d  decreases. 

SWAP  was  used  to  calculate  K  on  the  spall 
plane  for  overdriven  spall  for  shots  13,  24,  63,  and  223  (Reference 
19).  Values  of  K  were  computed  up  to  the  genesis  of  the  spall 
signal.  These  values  were  consistently  about  one-half 
agreement  with  the  STEALTH  calculation.  However,  the  discrepancy 
in  K  values  for  threshold  and  overdriven  spall  is  probably  of 
little  practical  importance.  The  delay  in  the  genesis  of  the  spall 
signal,  if  spall  is  assumed  to  occur  at  K  =  is  only  about 

20  ns.  This  is  comparable  to  the  uncertainty  in  wave  arrival 
times  that  is  caused  by  the  assumption  of  perfectly  elastic/plastic 
behavior  embedded  in  SWAP. 

While  the  Tu ler-Butcher  model  can  be  used  to 

predict  spall,  it  is  probably  not  physically  correct.  The  spall 

data  indicate  that  void  nucleation  begins  at  a  =17  Kbar,  not  10 

s 

Kbar  in  the  plate  impact  tests  (see  Section  2,3.3  a).  The  spall 

signals  show  the  a  increases  with  impact  velocity. 

s 

b.  OFHC  copper 

Bless  et  al. ,  (Reference  6)  summarized  the 
available  spall  results  on  fully  annealed  copper.  A  STEALTH 
simulation  of  their  test  No.  560  was  carried  out  (4mm  flyer,  6mm 
target,  impact  velocity  209  m/s).  The  spall  parameters  for  the  two 
criteria  were  obtained  from  the  simulation.  The  simulations  and 
VISAR  data  are  compared  in  Figure  106.  The  corresponding  values 
are  given  in  Table  10.  The  value  was  assumed  to  be  egual  to  the 
Hugoniot  stress  3.6  Kbar,  and  A  was  arbitrarily  assigned  a  value 
of  2.  The  value  for  K^  in  equation  49  was  obtained  by  matching 
the  arrival  times  of  spall  signals  in  the  simulation  and  test. 

It  was  found  that  the  spall  stress  computed 
from  equation  13  was  always  between  11  and  12  Kbar.  A  much 
higher  value  of  spall  strength,  31  Kbar,  was  required  to  obtain 
agreement  between  velocity  profiles  calculated  with  SWAP  and 
observed  with  the  VISAR  for  shots  well  above  the  spall  threshold. 


VELOCITY  (KM/SEC) 


TIME  (MICROSECONDS) 


Figure  106.  Spall  Simulation  -for  OFHC  Copper  Target  and 
1020  Flyer  from  STEALTH  Using  Tuler-Butcher 
Spall  Model  as  Compared  with  VISAR  Data. 


191 


TABLE  10 


CRITICAL  SPALL  STRESS  AND 
TULER-BUTCHER  MODEL  PARAMETERS 


TULER-BUTCHER  MODEL 

SIMPLE 

SPALL 

MODEL 

ACTUAL 

SPALL 

STRESS 

Material 

a 

o 

(Kbar) 

X 

K, 

Xc 

a 

c 

(Kbar) 

(Kbar) 

1020 

STEEL 

10 

2 

12.5 

30 

17 

OFHC 

3.6 

2 

10.0 

24.7 

13 

C1008 

14.0 

2 

0.38 

20.0 

18 

HYIOO 

15.7 

2 

61.0 

56.8 

27 

7039-T64 

8.6 

2 

3.0 

18.2 

13 

A  constant  spall  stress  model  in  SWAP  could  not  explain 
observations  near  the  spall  threshold  (Reference  46). 

Several  SWAP  calculations  were  conducted  to 
match  the  observed  variation  of  critical  impact  velocity  with 
flyer  plate  thickness  (Reference  6).  The  variation  of  the  spall 
threshold  with  flyer  plate  thickness  could  be  described  with  a 
Tuler-Butcher  model  with  1=2  and  =  7.5  Kbar. 

As  noted  previously,  the  spall  signal  in  1020 
steel  appeared  to  be  associated  with  void  nucleation.  If  that  is 
true  for  OFHC  copper  also,  and  if  the  spall  develops  gradually, 
then  the  interpretation  that  11  Kbar  is  the  stress  at  which 
damage  initiates  may  be  physically  correct.  However,  when  using 
a  computer  code  in  which  complete  rupture  occurs  instantaneously, 
accurate  predictions  of  threshold  are  obtained  with  the  derived 
Tuler-Butcher  parameters,  and  predictions  of  spall  signals  for 
high  impact  velocities  are  obtained  with  a  spall  stress  of  31 
Kbar.  These  results  are  discussed  more  extensively  by  Bless  and 
Paisley  (Reference  46). 

An  error  in  the  SWAP-based  Tuler-Butcher 
predictions  was  revealed  during  experiments  with  double  flyer 
plates  (described  in  Section  2.4.1).  In  those  experiments, 
aluminum  flyer  plates  and  copper  targets  were  employed.  SWAP 
simulations  were  used  to  predict  spall  in  the  copper.  However, 
it  was  found  that  the  Tuler-Butcher  parameters  given  above 
substantially  underestimated  the  impact  velocity  required  to 
cause  damage.  The  cause  of  this  discrepancy  has  not  yet  been 
resolved.  It  may  be  related  to  the  strain  rate  on  the  spall 
plane.  The  strain  rate  varies  when  a  given  peak  stress  value  is 
obtained  by  changing  flyer  plate  materials  and  velocities. 

Strain  rate  does  not  explicity  appear  in  the  Tuler-Butcher 
formulation.  Therefore,  the  spall  criterion  for  OFHC  copper  must 
be  used  with  caution  for  situations  that  differ  greatly  from 
those  in  the  copper/copper  plate  impact  experiments. 
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c. 


1008  steel 


The  3min  thick  flyer  plates  induced  spall 

(defined  as  50  percent  voids  in  the  spall  plane)  at  a  critical 

velocity  of  168  iti/s.  This  corresponds  to  a  stress  of  31  Kbar. 

However,  a  ,  was  found  to  be  18  Kbar.  A  STEALTH  simulation  of 
s 

shot  678  reproduced  the  observed  spall  signal  with  a  simple  stress 
model  with  =  20  Kbar.  These  observations  suggest  that  in  these 
three  shots,  spall  developed  over  a  finite  period  of  time.  The 
spall  signal  is  associated  with  void  nucleation,  and  this  takes 
place  at  about  20  Kbar. 

A  STEALTH  simulation  was  also  made  with  the 
Tuler-Butcher  model  for  test  678.  The  exponent  X  was  assumed  to 
be  *2',  as  for  1020  steel.  The  stress  a  was  set  equal  to  o  , 

O  HhiL 

14  Kbar.  The  parameter  K  was  adjusted  in  the  simulations  until 

A 

the  experimental  results  were  matched.  The  calculated  free 
surface  velocity  of  the  target  is  compared  with  the  VISAR  data  in 
Figure  107.  The  corresponding  stress  history  at  the  spall  plane 
is  also  shown.  When  K  =  K  ,  spall  occurred  and  the  stress  in 

A  AC 

the  spall  plane  was  instantaneously  set  to  zero  in  the 

simulation.  The  free  surface  velocity  history  and  spall  plane 

history  from  these  two  models  are  almost  identical  as  can  be  seen 

from  Figure  107.  Neither  model  precisely  predicts  the  post  spall 

behavior,  but  they  are  probably  accurate  enough  for  most 

engineering  applications.  Since  the  simple  spall  criterion  is 

easier  to  use  than  the  Tuler-Butcher  model,  it  is  recommended. 

However,  it  must  be  remembered  that  these  models  do  not  describe 

damage  development.  Their  utility  is  for  prediction  of  the  onset 

2 

of  tensile  failure.  The  value  of  K^^  was  only  0.38  Kbar  ys. 

Thus,  the  material  spalled  very  near  o^.  Probably,  the  Tuler- 
Butcher  model  would  have  worked  equally  well  with  a  =  a  =18 
Kbar,  since  the  stress  on  the  spall  plane  falls  extremely  rapidly 
when  the  release  waves  arrive. 

The  Tuler-Butcher  fracture  parameters 
estimated  from  simulations  of  shot  678  were  used  to  simulate  test 
#687  which  was  at  a  lower  impact  velocity  (u^  =  134  m/s).  The 
VISAR  traces  for  shot  687  indicated  a  "spall  signal".  The 
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STEALTH  simulation  reasonably  reproduced  these  features,  as  shown 

in  Figure  108.  Metallographic  and  SEM  examination  showed  the 

damage  in  shots  685  and  687  (same  velocity)  consisted  of  isolated 

voids,  rather  than  complete  spall.  In  summary,  based  on  these 

limited  data  (only  one  thickness  flyer  plate),  spall  initiation 

can  be  described  by  a  simple  spall  stress  ,  o  =  20  Kbar,  or  a 

Tuler-Butcher  model  with  A  =  2,  a  =15  Kbar,  and  K  =0.38 

2  o  crit 

Kbar"^  ys.  Complete  rupture  on  the  spall  plane  requires  a  stress 

of  31  Kbar  if  damage  development  is  not  simulated.  However,  the 

value  of  the  stress  actually  associated  with  void  formation  is 

probably  18  Kbar. 

d.  HYIOO  steel 


A  STEALTH  simulation  of  test  #680  on  HYIOO 

was  carried  out  with  the  two  spall  criteria.  The  overlapping  results 

are  shown  in  Figure  109.  The  critical  spall  stresses  and  the 

Tuler-Butcher  fracture  parameters  that  were  assumed  in  these 

simulations  are  given  in  Table  10.  The  a  value  was  assumed  to  be 

o 

equal  to  the  Hugoniot  stress,  o  obtained  from  the  plate-impact 

nhilj 

tests.  (Preliminary  values  were  used  that  differ  slightly  from 
those  of  Table  2.  )  The  A  was  arbitrarily  taken  to  be  2.  The 
values  of  the  critical  damage  parameter,  K^^  are  adjusted  to 
predict  the  experimental  results. 

e.  7039-T64  aluminum 


The  test  #682  on  7039-T64  was  simulated 
through  STEALTH  with  a  simple  material  model.  The  elastic-plastic 
models  with  Y  =4.6  for  the  flyer  (1020  steel)  and  Y  =4.2  for 
the  target  (7039-T64),  were  used  in  the  simulation.  Since  the 
release  characteristics  depend  only  on  the  Y^,  the  VISAR  and  the 
simulation  compare  extremely  well  in  Figure  110.  Two  spall 
criteria  were  used  in  the  model  calibration.  The  critical  spall 


stress  ( a  )  and  the  Tuler-Butcher  fracture  parameters  that  were 
s 

assumed  in  the  simulations  are  given  in  Table  10.  The  o  value 

^  o 

was  assumed  to  be  equal  to  For  aluminum,  A  was  also 

rltiJLt 

arbitrarily  taken  to  be  2.  The  value  of  the  critical  damage 
parameter,  K^^  was  adjusted  to  predict  the  experimental  results. 
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Though  the  VISAR  data  are  well  matched  in 
these  predictions,  the  estimated  fracture  model  parameters  are 
probably  not  unique  and  may  lack  generality.  (Experiments  would 
need  to  be  conducted  with  other  flyer  plate  thicknesses  to  check 
this  point. )  In  these  models,  failure  is  assumed  to  occur  almost 
instantaneously.  Models  which  predict  failure  successfully  under 
high  impact  velocity  conditions  often  cannot  predict  failure 
under  threshold  conditions  (Reference  19).  Zukas  (Reference  1) 
has  documented  well  the  varying  nature  of  several  other  failure 
models  and  their  inherent  inadequacies  in  supporting  numerical 
investigations  of  impact  events. 


SECTION  4 
SUMMARY 


High  strain-rate  material  characterization  was  carried 
out  for  1020,  C1008,  and  HYlOO  steels,  OFHC  copper,  7039-T64 
aluminum,  and  BeO  ceramic.  Tests  were  conducted  under  (1) 
uniaxial  stress  (compression  and  tension);  (2)  triaxial  stress 
(necking  specimen);  and  (3)  one  dimensional  strain  conditions. 
Split  Hopkinson  bar  tensile  and  compression  tests  were  used  to 
achieve  uniaxial  stress  states.  High  speed  photography  of  necking 
SHB  tensile  test  specimen  provided  data  under  triaxial  stress 
states  and  large  strain  conditions.  Plate  impact  tests  were  used 
to  achieve  one  dimensional  strain  conditions.  Table  11  summarizes 
the  principal  diagnostics,  conditions  and  the  experimental 
parameters  involved  in  these  three  test  configurations. 

The  wide  range  of  experimental  data  obtained  from  the 

various  tests  were  used  for  high  strain-rate  material  modeling. 

The  state-variable-based  visco-plastic  theory  of  Bodner  and  Partom 

was  used  to  describe  material  behavior  under  a  broad  range  of 

stress,  strain,  and  strain-rate  conditions.  Computer  programs  and 

additional  subroutines  for  the  STEALTH  finite  difference  code  were 

developed  to  evaluate  material  model  parameters.  A  unique 

iterative  algorithm  was  formulated  to  use  both  SHB  and  plate 

impact  test  results.  The  parameters  for  the  B-P  model,  D  ,  n,  m  , 

o  o 

m, ,  a  ,  Z  ,  and  Z  were  obtained  for  all  five  metals  tested  with  D 
loo  o 

assumed  constant.  The  three  steels  and  aluminum  were 
characterized  using  only  ,  n,  m^,  Z^ ,  and  Z^.  However  OFHC 
copper,  due  to  its  strong  strain-hardening  behavior,  was  modeled 
with  additional  constants  m^  and  a. 

Dynamic  tensile  failure  in  the  selected  metals  was 
characterized.  Since  the  objective  was  to  provide  the  capability 
CO  predict  failure  in  engineering  calculations,  simple  failure 
models  were  used.  Spall  failure  was  characterized  with  1)  a 
simple  time  independent  critical  spall  stress  criterion,  and  2) 
the  Tuler-Butcher  time  dependent  damage  model.  The  spall  models 


fO 

<D 

C 

•H 

e 

Sh 

(U 

4J 

Q) 

Q 

g 

0 

0) 

•H 

U 

+j 

<0 

Cu 

o 

+> 

(l) 

g 

0) 

+J 

& 

g 

u 

1 — ! 

g 

(C 

m 

+> 

rH 

0 

u 

kk 

M-l 

iH 

U 

p 

g 

o 

(0 

-P 

e 

g 

• — 

H 

d 

g 

fit 

O 

G 

G 

>H 

M 

CO 

P3 

w 

w 

D 

a 

o 

H 

fH 

CO 

1 

o 

X 

CO 

(H 

g 

u 

— ‘ 

0 

M 

o 

p 

EH 

<u 

CO 

-P 

ip 

rH 

s 

o 

1 

EH  2 

(Ci 

iH 

Ai 

Td 

g 

g 

2  K 

X 

CD 

CD 

+j 

W  O 

o 

CD 

-p 

VO 

g 

VO 

S  O 

g 

-P 

g 

(d 

o 

p 

o 

M  X 

•H 

p 

iH 

g 

1 — 1 

X  Cl^ 

(d 

>1 

g 

e 

rH 

rH 

D 

t 

'H 

Ip 

1 

rH 

1 

r— i 

W  Ui 

-P 

(d 

CD 

g 

|<  H 

CO 

CD 

■«r 

o 

c» 

0 

in 

w 

W  X 

+J 

o 

o 

o 

g 

o 

hi 

a  eh 

CO 

iH 

2 

rH 

p 

rH 

m 

m  x 

Eh 

&H  H 

^  Q 

g 

W 

O 

X  W 

1 — 1 

■H 

H  X 

m 

■p 

w  g 

X 

c 

g 

g  o 

Eh 

o 

0 

P  'H 

g 

cn 

'H 

43 

g 

d  g 

O 

g 

■p 

(d 

-H 

g  g 

■iH 

'H 

K 

cd 

g 

0  g 

g 

g 

a 

g 

0 

p 

U  -P 

g 

p 

H 

p 

+J 

HJ 

g 

p 

X 

0 

t) 

G 

CD 

-p  g 

p 

g 

m 

>1 

g  g 

& 

Q) 

Q 

P 

Q 

P  43 

g 

P 

Q 

1 

-g 

1 

1 

■H  -P 

o 

1 

r-H 

G 

iH 

rH 

d 

o 

I — 1 

g 

o 

* 

-H 

4J 

CD 

■p 

CJ 

g 

g 

CD 

u 

g  g 

+j 

CD 

0 

CD 

g 

d  p 

g 

g 

"H 

P 

d 

g  g 

1 — 1 

■p 

CO 

d 

g 

H  g 

d* 

(d 

g 

g 

H 

g 

+j 

u 

(D 

O 

P  o 

g  o 

(d 

Eh 

u 

CD 

0 

H  g 

a 

1 

■P 

H  CO 

43  d 

d 

« 

ffl 

g 

>1  1 

g  g 

2 

ffi 

X 

rH 

g  o 

O  H 

CO 

CO 

d 

Eh  • — 

p 

202 


Technique  that  was  developed,  but  has  not  yet  been  used  for  evaluating  material  parameters. 


were  used  only  with  a  rate-independent,  elast ic-per fectly  plastic 
material  description. 

The  following  sections  summarize  the  critical  test  data  and 
the  model  parameters  in  the  high  strain-rate  material  character¬ 
ization  of  1020,  C1008,  HYIOO  steels,  OFHC  copper,  7039-T64 
aluminum,  and  BeO  ceramic. 

4.1  SUMMARY  OF  CRITICAL  TEST  DATA 

The  critical  test  data  from  the  various  experiments  are 
presented  in  terms  of  figures  and  tables  for  all  six  materials. 
Based  on  these  results,  certain  observations  on  rate  dependency, 
strain  at  failure,  stress  at  failure,  and  spall  characteristics 
were  made  and  are  summarized  in  the  following  paragraphs. 

a .  1020  steel 

This  steel  showed  considerable  rate  dependency 

-4 

between  the  strain  rates  of  10  /s  and  400/s  as  can  be  seen  from 

Figure  111.  Since  the  strain  rates  associated  with  the  HEL  were 

not  experimentally  obtained,  the  HEL  is  represented  by  a  scatter 

4  5-1 

band  between  the  strain  rates  of  10  and  10  s  in  Figure  111. 
Between  400  and  1200/s,  there  was  less  rate  dependency.  However, 
tjrends  could  not  be  quantified  due  to  test  to  test  variability. 
Compression  data  at  2000  s  ^  agreed  with  tensile  data  at  1200  s  ^ 
and  HEL  data,  which  implies  rate  independency  and  isotropy.  Key 
data  are  displayed  in  Figure  5. 

The  effective  stress  vs.  effective  plastic  strain 
data  obtained  from  necking  specimens  are  shown  in  Figure  41.  The 
increase  in  stress  due  to  an  increase  in  strain-rate  during  the 
necking,  is  not  significant  for  this  material.  The  stress  at 
large  strains  and  at  higher  strain-rates  agreed  with  the  stress  at 
iOOO/s,  indicating  that  this  material  is  relatively  rate- 
independent  above  a  strain  rate  of  approximately  400/s. 

The  maximum  stress  obtained  at  the  end  of  the  first 
pulse  from  the  SHR  test,  along  with  the  strains  at  onset  of 
necking  and  failure  are  summarized  in  Table  3. 
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gure  111.  Strain  Rate  vs.  Stress  Behavior  for  1020  Steel. 
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b. 


OFHC  copper 


It  is  well  known  that  OFHC  copper  is  a  highly  rate 
dependent  material.  SHB  test  data  clearly  showed  this  behavior. 

In  Figure  112  stress  vs.  strain-irate  at  various  strain  levels  is 
plotted.  The  quasi-static  flow  stresses  are  also  included.  This 
material  exhibited  a  strong  strain-hardening  even  at  large  strain 
values.  The  strain  hardening,  ,  is  an  increasing  function  of 
strain  rate  over  the  range  of  strain  and  strain  rate  tested. 
Comparison  of  tension  and  compression  data  shows  this  copper  is 
anisotropic. 

The  results  of  various  SHB  tensile  tests  are  shown 
in  Figure  11.  Necking  did  not  occur  during  the  first  tensile 
pulse,  so  the  effective  stress  at  large  local  strains  could  not  be 
obtained.  The  failure  parameters  obtained  from  the  photographs 
are  summarized  in  Table  3.  Various  HEL  and  spall  parameters  from 
earlier  work  are  summarized  in  Table  12. 

c.  C1008  steel 

Quasi-static  and  SHB  tensile  and  compression  tests 
were  conducted.  Stress  vs.  strain-rate  plots  at  different  strain 
rate  levels  are  shown  in  Figure  15.  Based  on  these  data,  it  is 
concluded  that  C1008  steel  is  a  rate  dependent  material.  When 
the  initial  yield  stresses  at  various  intermediate  strain-rate 
levels  were  compared  with  the  (at  very  high  strain-rate 

level),  the  rate  dependency  was  quite  evident.  Figure  113 
summarizes  the  various  stress-strain  data  obtained  from  both  the 
quasi-static  and  SHB  tests. 

Necking  did  not  occur  during  the  first  load  cycle. 
Hutchinson  and  Neal  (Reference  50)  proposed  that  rate  dependency 
delays  neck  formation.  This  theory  is  consistent  with  our  obser¬ 
vation  on  OFHC  copper  and  C1008  steel.  Failure  strains  were 
calculated  from  post  test  measurements  and  are  reported  in  Table  3. 

Several  plate  impact  tests  were  conducted  and  free 
surface  data  were  obtained.  Rounding  and  ramping  were  observed 
in  most  free  surface  velocity  records.  Code  simulations  showed 


TABLE  12 


OBSERVED  HEL  AND  SPALL  THRESHOLDS  FOR  CRACK  FORMATION 


Material 

** 

HEL  * 

( Kbar ) 

Spall  ** 

Impact  Velocity 
(M/S) 

Spall 

Stress,  a 
(Kbar)  ^ 

C1008 

11.3 

168  +  10 

18 

HYIOO 

15.5 

190  +  10 

27 

7039-T64 

7.7 

>162 

13 

1  Khar  =  14504  psi 

=  100  MPA 

9  2 

=  10  dynes/cm 

*  * 

3iPin  flyer  plates 
Gmin  target  plates 


EFFECTIVE  STRESS  CKBAR) 


* 


Figure  112 


Strain  Rate  vs.  Stress  Behavior  for  OFHC  Copper. 


E  STRESS  CKBAR) 


that  this  behavior  is  caused  by  strong  dependence  of  flow  stress 
on  strain  rate.  (This  effect  is  masked  in  copper  by  large  strain 
hardening. )  In  Table  12,  results  from  the  plate  impact  tests  are 
summarized  in  terms  of  values  for  stresses  and  velocities  at  HEL 
and  spall  threshold. 

d.  HYIOO  steel 

The  experiments  on  HYIOO  steel  indicate  that  it 
is  moderately  strain-rate  sensitive.  Quasi-static  and  several  SHE 
test  results  are  summarized  in  stress  vs.  strain-rate  plots  at 
various  strain  levels  in  Figure  114.  The  Y^  obtained  from  the 
plate  impact  test  results  is  also  included  in  this  figure.  The 
increase  in  the  initial  flow  stress  due  to  strain-rate  is 
moderate  for  this  material.  The  stress-strain  data  (average 
curve  fits  of  the  test  data)  from  the  SHE  tensile  tests  and  the 
guasi-state  test  are  summarized  in  Figure  19.  The  slight  rate 
dependency  of  this  material  can  be  seen  from  this  figure. 

Three  necking  experiments  were  conducted  on  HYIOO 
steel.  The  effective  stress  vs.  effective  plastic  strain  deduced 
from  the  photographs  using  Eridgman  analysis  is  plotted  in  Figure 
46.  The  value  of  effective  stress  at  60  percent  strain  was 
around  12  Kbar.  Comprehensive  failure  parameters  obtained  from 
these  necking  tests  are  summarized  in  Table  3.  Several  plate 
impact  tests  were  conducted.  Excellent  VISAR  data  were  obtained 
and  the  results  are  summarized  in  Table  12. 

e .  7039-T64  aluminum 

The  results  in  terms  of  stress  vs.  strain-rate  from 
the  SHE  and  quasi-static  tests  are  plotted  in  Figure  115.  Y^ 
derived  from  the  plate  impact  tests  is  also  included  in  the 
figure.  The  tests  clearly  indicate  that  7039-T64  aluminum  is  a 
rate  insensitive  material.  Sample  stress-strain  data  are 
displayed  in  Figure  22. 

Three  necking  experiments  were  conducted.  The 
aluminum  samples  did  not  begin  to  neck  during  the  first  pulse. 
During  the  second  pulse,  the  necking  began  at  a  strain  of  25 
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Strain  Rate  vs.  Stress  Behavior  of  7039-T64 
Aluminum. 


percent  and  the  sample  failed  rapidly  with  little  additional 
necking.  The  failure  strain  listed  in  Table  3  was  obtained  from 
post-test  measurements.  Reshocks  in  this  material  traveled  at 
the  elastic  wavespeed.  The  plate  impact  test  results  are 
summarized  in  Table  12. 

f .  BeO  ceramics 

The  Hugoniot  of  polycrystalline  BeO  with  0.6  to  1.6 
percent  porosity  was  measured  to  110  Kbar.  The  data  are  presented 
in  Figure  64.  The  measured  HEL  values  were  extrapolated  to  96  Kbar 
for  zero-porosity  ceramic.  The  spall  stress  for  elastic  shocked 
material  was  about  1.8  Kbar.  The  spall  stress  of  material  shocked 
above  the  HEL  was  negligible.  In  material  experiencing  only  elastic 
shocks,  spall  seems  to  occur  over  a  region,  rather  than  on  a  well 
defined  plane  (Figure  67). 

4.2  SUMMARY  OF  MODEL  PARAMETERS 

In  this  section,  the  results  from  the  Bodner-Partom 

material  modeling  are  summarized  for  the  five  metals.  The 

calibrated  model  parameters  are  described  in  Table  10.  The 

parameter  n  indicates  the  level  of  strain-rate  sensitivity.  Its 

value  is  inversely  proportional  to  the  level  of  rate  dependency. 

For  example,  a  value  of  0.3  indicates  significant  rate 

sensitivity  and  a  va]ue  of  4  indicates  negligible  sensitivity. 

The  parameter  m  (or  m^)  indicates  the  level  of  strain  hardening. 

A  higher  value  for  m  indicates  strong  strain  hardening.  The 

difference  in  the  values  of  Z  and  Z,  indicates  the  difference 

o  1 

between  initial  yield  stress  and  stress  at  large  strain. 

For  each  material,  the  model  simulations  and  the  corre¬ 
sponding  test  results  are  summarized  in  various  figures.  Two 
spall  criteria,  (critical  spall  stress  and  Tuler-Butcher  model) 
were  used  in  the  spall  modeling.  The  model  parameters  were 
obtained  from  STEALTH  simulations.  The  results  are  summarized  in 
Table  10. 
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1020  steel 


a . 

Based  on  the  experimental  results,  this  material  is 
characterized  as  rate  independent  (for  the  strain  rates  of 
interest  here) .  The  interpretation  is  based  on  SHB  and  plate 
impact  data  for  1020  steel  targets  and  for  C1008,  HYIOO,  and 
7039-T64  aluminum  targets  struck  by  1020  steel  flyer  plates.  To 
accurately  simulate  release  characteristics  with  STEALTH,  it  was 
necessary  to  describe  1020  steel  with  =  4.6  Kbar,  the  value 
obtained  in  SHB  tests.  When  Bless  et  al.,  (Reference  44)  modeled 
the  behavior  of  1020  steel  with  the  HULL  code,  they  had  to  assume 
a  stress-strain  model  which  is  consistent  with  the  present  SHB 
test  results,  providing  additional  evidence  that  1020  steel  is 
relatively  rate  independent.  The  Bodner-Partom  material  model 
parameters  are  given  in  Table  9.  This  table  also  contains  the 
values  of  elastic  constants  obtained  from  the  literature.  The 
program  BPSOLVE  was  used  to  simulate  the  tensile  and  compression 
SHB  tests  with  the  BP  parameters  (Table  9).  The  agreement  with 
the  sets  of  SHB  data  was  excellent.  Representative  comparisons 
are  shown  in  Figure  83.  STEALTH  was  used  with  the  same 
parameters  to  simulate  plate  impact  tests  (Figure  85). 

b.  OFHC  copper 

This  rate  sensitive  material  could  not  be  modeled 

with  only  Z  ,  Z,,  n,  m  ,  and  D  .  Additional  constants  m,  and  a 
o  1  o  o  1 

were  needed  due  to  its  strong  strain  hardening.  Large  differences 

between  the  initial  yield  stress  and  the  flow  stresses  at  high 

strains  were  observed.  The  values  of  Z  (=9)  and  Z,(=72)  clearly 

o  i 

indicate  this  behavior.  The  model  constants  are  given  in  Table  9. 
The  model  predictions  compared  extremely  well  with  the  test  results 
as  shown  in  Figures  86  and  88.  The  spall  model  parameters  obtained 
through  SWAP  simulations  are  summarized  in  Section  3.4. 

c .  C1008  steel 

The  ability  of  the  BP  model  to  describe  the  SHB  and 
plate  impact  data  was  fully  satisfactory.  The  parameters  are 
summarized  in  Table  9.  The  relatively  small  value  of  n  (0.4) 
indicates  significant  rate  sensitivity.  Use  of  the  simple  spall 
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model  was  adequate  to  simulate  the  observed  free  surface  velocity 
histories.  However,  data  are  available  only  for  one  flyer  plate 
thickness,  so  the  generality  of  the  spall  criteria  could  not  be 
fully  assessed.  This  is  also  true  for  HYIOO  and  7039-T64 
aluminum. 

d.  HYIOO  steel 

From  the  experimental  results,  HYIOO  steel  was 
considered  to  be  a  moderately  rate-sensitive  material.  This  is 
reflected  by  the  value  of  the  BP-model  parameter  n  (1.2)  as  given 
in  Table  9.  Using  the  model  parameters  in  Table  9  for  HYIOO 
steel,  simulations  of  SHB  tests  and  the  plate  impact  tests  were 
successfully  made  and  the  results  are  summarized  in  Figures  100 
and  102.  The  spall  parameters  for  the  two  spall  criteria  are 
summarized  in  Table  10. 

e .  7039-T64  aluminum 

The  Bodner-Partom  model  constants  were  obtained  from 
the  SHB  and  plate  impact  test  data.  Since  this  material  is  a  rate 
independent  material,  the  model  parameter  n(=4)  takes  a  large 
value.  A  convenient  feature  of  Bodner-Partom  visco-plastic  model 
is  its  ability  to  describe  even  a  rate  independent  solid  through 
suitably  selected  material  constants.  The  algorithm  developed  to 
describe  BP-model  in  STEALTH  is  then  valid  for  materials  with  and 
without  rate-sensitivity.  Similarly,  both  elastic-plastic  strain 
hardening  and  elastic-perfectly  plastic  material  can  be  described 
by  the  same  algorithm.  The  BP-model  parameters  enabled  successful 
simulations  of  both  SHB  tests  and  plate  impact  tests. 

4.3  RECOMMENDATIONS 

The  following  recommendations  on  experimental  techniques  and 
modeling  are  made  based  on  the  results  presented  in  this  report. 

4.3.1  Experimental  Technique 

The  expanding  ring  test  technique  should  be  refined 
and  applied  to  ductile  metals.  Data  from  this  test  should  be  used 
to  validate  the  Bodner-Partom  parameters  for  very  high  tensile 


strain  rates,  larger  strains,  and  varying  strain  rates.  This  test 
has  been  developed  in  previous  work  by  Warnes  (Reference  51), 

Bless  et  al. ,  (Reference  6)  and  Rajendran  and  Fyfe  (Reference  52). 
It  should  not  be  difficult  to  develop  procedures  to  routinely 
obtain  reliable  data  with  which  constitutive  models  can  be 
evaluated . 

The  double  flyer  plate  technique  should  be  applied 
systematically  to  metals  and  ceramics.  It  should  be  used  to 
obtain  data  for  the  strength  and  modulus  of  damaged  material  in 
or  near  the  spall  plane.  In  order  to  be  most  useful,  the  double 
flyer  technique  should  be  combined  with  soft  recovery  techniques. 

The  pressure-shear  impact  technique  developed  by 
Kim  and  Clifton  (Reference  53)  should  be  used  to  obtain  data  for 
validating  material  modeling  under  combined  normal  and  shear 
stress  state  conditions.  Experiments  in  which  conventional  spall 
failure  occurs  will  provide  data  for  the  strength  of  spalled 
material.  Experiments  in  which  failure  occurs  behind  the  shear 
wave  will  show  how  failure  depends  on  shear  stress.  These 
experiments  should  be  conducted  for  values  of  a  cc/ ^  that  are 
typical  in  impacts  of  engineering  or  military  importance. 

4.3.2  Material  Modeling 

The  evaluation  of  constitutive  model  parameters  must 
include  the  data  obtained  from  the  varieties  of  tests  so  that  the 
calibrated  model  could  predict  the  flow  stresses  under  different 
extremes  of  stress,  strain  and  strain-rate  states. 

More  work  is  needed  to  develop  a  comprehensive 
failure  criterion  that  is  of  sufficient  generality.  The  spall 
criteria  used  in  the  present  report,  for  example,  will  not  in 
general  describe  failure  in  SHB  tests  and  vice  versa.  Continuum 
mechanics  based  models  that  incorporate  fundamental  micro¬ 
mechanisms,  should  be  exercised  in  finite  difference  codes  for 
meaningful  predictions. 

The  failure  models  that  are  currently  available  in 
literature,  must  be  tested  for  their  generality.  Yield  models  of 


damaged  materials  which  explicitly  include  porosity,  the  second 
deviatoric  stress  invariant  and  the  mean  stress,  must  be  validated 
through  data  obtained  from  the  double  flyer  plate  test.  The  spall 
data  from  the  pressure-shear  experiment  will  be  useful  to  model 
the  influence  of  shear  stresses  on  damage  nucleation,  growth,  and 
failure  models. 

Additional  plate  impact  experiments  with  thicker 
flyer  plates  (o/lOmm)  should  be  conducted  on  1020,  C1008,  and  HYlOO 
materials.  Data  should  be  obtained  for  precursor  decay.  The 
anomalous  ramping  in  C1008  steel  should  be  verified  and 
investigated. 

Finally,  the  developed  techniques  must  be  used  to 
characterize  other  materials  of  engineering  importance,  for 
example;  high  density  penetrator  alloys,  structural  metals, 
metals  used  in  explosive  forming,  composite  materials,  and 
ceramics  such  as,  BeO  and  TiB2. 
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Flow  Diagram  for  Bodner-Partom  Model  in  STEALTH. 
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Flow  Chart  for  the  Subroutine  SIJNEW, 


A- 4 


•f 


'i 


SDOTtJ  UH*  ch*  (S*")^ 

CO  eoapuco  now  dovtocorte 
icrooo  grodUoc  ••claacoo. 


Ij 


« 


A4.  Flow  Chart  for  the  Subroutine  SIJEST.  Improved 
Sij  Values  are  Estimated  in  this  Routine. 
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A5.  Flow  Chart  for  the  Subroutine  SDOTIJ.  Sij  are 
Computed  in  this  Routine. 
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Flow  Chart  for  the  Subroutine  EPIJ.  e?  . 
Computed  in  this  Routine. 


A- 7 


A7. 


r 


L 


% 


> 


Flow  Chart  for  the  Subroutine  SIJDIF.  Convergence 
Check  of  the  Iterative  Scheme. 
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